Abstract 



We generalize the Umbral Calculus of G-C. Rota [40] by studying not only sequences of polynomials and 
inverse power series, or even the logarithms studied in [24], but instead we study sequences of formal expres- 
sions involving the iterated logarithms and x to an arbitrary real power. 

Using a theory of formal power series with real exponents, and a more general definition of factorial, 
binomial coefficient, and Stirling numbers to all the real numbers, we define the Iterated Logarithmic Algebra 
X. Its elements are the formal representations of the asymptotic expansions of a large class of real functions, 
and we define the harmonic logarithm basis of X which will be interpreted as a generalization of the powers 
x n since it behaves nicely with respect to the derivative. 

We classify all operators over X which commute with the derivative (classically these are known as shift- 
invariant operators), and formulate several equivalent definitions of a sequence of binomial type. We then 
derive many formulas useful towards the calculation of these sequences including the Recurrence Formula, 
the Transfer Formula, and the Lagrange Inversion Formula. Finally, we study Sheffer sequences, and give 
many examples. 



L'algebre des Logarithmes Iteres 

On generalise ici le calcul ombral de G-C Rota en etudiant non seulement les suites de polynomes, de 
series a exposants entier negatif, ou encore de logarithmes [24], mais d'une fa con plus generale on s'interesse 
aux suites formelles des logarithms iteres de x et de x a un exposant reel quelconque. 

On s'appuie sur la theorie des series de puissances formelles a exposants reels, sur une definition generale 
des factorielles, des coefficients du binomes, et des nombres de Stirling a tous les nombres reel pour definir 
l'Algebre de Logarithmes Iteres X. Ccttc algebre a comme elements les representations asymptotiques de 
beaucoup de fonctions reelles par rapport de Pechelle des monomes dans les Logarithmes Iteres. On definit 
alors une base de X de logarithms harmoniques, lesquels pourront etre considered comme generalisant les 
puissances x n , puisqu'ils se comportent comme x n lors de la derivation. 

On caracterise tous les operations sur X qui commutent avec l'operateur de derivation (ce sont usuellement 
les "operateurs invariants par translation" ) , ce qui nous permet de formuler plusieurs definition equivalentes 
des suite de series logarithmiques de type binomial. On en deduit de nombreuses formules utiles au calcul 
de ces series, en particulier la formule de recurrence, la formule de transfert, et la formule d'inversion de 
Lagrange. Enfin, on etudie les suites du Sheffer, ainsi que de nombreux exemples. 
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Chapter 1 



Introduction 



Over the years, mathematicians have studied many special sequences of functions especially polynomials. 
These sequences were found to have many important similarities; however, it was not until G-C. Rota's 
Umbral Calculus [40] that this notion was formalized. Only could one study under a single theory most 
of the important sequences of polynomials — for example, x n , Abel, lower factorial, upper factorial. These 
sequences are all sequences of binomial type. 

Nevertheless, yet more work needed to be done, for the theory only applied to polynomials. For example, 
the theory of "factor sequences" was developed specifically to handle inverse formal power series. However, 
this theory was completely separate from the Umbral calculus of polynomials. An early version of this article 
was published under the title Formal Power Series of Logarithmic Type [24] . Not only did this paper allow 
us to simultaneously consider polynomials and inverse formal power series, but it also allowed us to consider 
expressions involving the logarithm logx. 

This article differs from [24] principally in that it considers two further generalizations to the logarithmic 
algebra C. By introducing the Iterated logarithms log logx, log log logx, and so on to the Logarithmic 
algebra, we have the discrete theory. Then by including x to any real power, we derive the continuous theory. 
Thus, the logarithmic algebra C is the set of all complex functions on the real numbers with asymptotic 
expansions in a neighborhood of +oo with respect to the ladder of comparison x a (logx) b • • •. With only a 
handful of obvious exceptions, all of the results of Formal Power Series of Logarithmic Type carry over into 
these new contexts with only minor changes. For example, we conclude that every sequence of polynomials 
binomial type can be uniquely extended to a pseudobasis of C called a Roman graded sequence (after Prof. 
Steve Roman). In fact, this sequence is usually easier to calculate than the sequence of binomial type itself 
using classical techniques. Thus, even when one is only interested in polynomials it still pays to introduce 
logarithms. 
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CHAPTER 1. INTRODUCTION 



1.1 Discrete and Continuous 



In this paper, all results and sections regarding the discrete iterated logarithmic algebra will be denoted 
"Discrete" and the results regarding the more general continuous iterated logarithmic algebra will be denoted 
"Continuous." Readers interested in only one of these theories may safely omit all material pertaining to 
the other. Sections and results numbered with a D and paragraphs starting with Discrete are relevant 
only to the discrete theory whereas sections and results numbered with a C and paragraphs starting with 
Continuous are relevant only to the continuous theory. Otherwise, the remainder of this article may be 
interpreted discretely by supposing that all variables a, b, c, . . . are integers, and that a, (3 are vectors of 
integers, or it may be interpreted continuously by supposing that all variables a,b,c,... are real numbers 
and that a, (i are vectors of real numbers. 

Digressions and all sections marked with the word "Appendix" or the letter "A" are independent of all 
later material, and are included for their own sake. 



Chapter 2 

D-Invariant Operators 



2.1 The Operator Topology 



It is a classical result that the algebra of formal differential operators X)«>o a «D™ acts on the vector space 
of polynomials. In view of the fact that the derivative D is invertible in X a for a ^ (0) (Proposition ??), 
we can define on I + the action of a more general class of differential operators, called Artinian operators. 
These are linear operators which commute with all powers of the derivative, and act on each level in a similar 
fashion; this notion will be made more precise later. 

The primary goal of this chapter is to classify all Artinian operators. We see (Theorem 2.4.1) that they 
are merely Artinian series in the derivative, and that the operator topology — which we are shortly to define — 
corresponds to the Artinian topology. This implies (see §2.4.4A) that all linear differential equations have a 
unique canonical solution. Finally, along the way we find the opportunity (in §2.3) to apply this theory to 
Real Analysis. 

We begin by defining a topology on the ring of continuous linear operators acting on formal power series 
of logarithmic type. 

Let J be a subspace of I. We say that a sequence {6 n )n>o of continuous linear operators of J into itself 
converges in the operator topology on J when for every p{x) e J the sequence (9 n p{x)) n >o converges in J . 
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Proposition 2.1.1 Let J = T, I + , orI a . Then the set of continuous linear operators in the operator 
topology of J is a complete topological K -algebra whose operations are given by: 



Proof: Let (6*„)„> and (<j) n ) n >o be convergent sequences of continuous linear operators on J. Thus, 
for any p(x) G J , the sequences (O n p(x)) n >o and (<fi n p(x))n>o are Cauchy. In particular, (0 n 4>kP(x)) n >o is 
Cauchy for any k > 0, and since Ok is continuous, (0k<p n p(x)) n >o is also Cauchy for any k > 0. Hence, 
{0 n <f>nP(x))n>o is Cauchy, and (0 n (f> n ) n > converges. 

((On + 4*n)p(x))n>v converges since J' is a topological space. Hence, (0 n + (j} n )n>o converges. 

The ring is complete since 6p(x) = lim„^ +00 n p(x) is the limit of a Cauchy sequence. E 

We shall write infinite series ^2 k>d Ok of operators, which are understood to denote the limits of their 
partial sums. 

We list below some notable operators and comment briefly on them: 



Example 2.1.1 (Derivative) The derivative D a is defined via Definition ?? on all of 1 when a is a 
nonnegative integer. Otherwise, it is defined on only the positive logarithmic algebra I + via 

Discrete Proposition ??. 

Continuous Definition ??. Note that the map a i— > D a is continuous in the operator topology. 



Example 2.1.2 (Shift Operator) For all complex numbers z, the shift operator, E z : 1 — > 1 is given by 
the sum 



(0<f>) P (x) 
(aO)p(x) 
(0 + 4>) P (x) 



0(<j>p(x)) 
a(0p(x)) 

(0p(x)) + (#(*)). 



D a XUx) = 



[b]l 



[b-a]\ 




n>0 



We see that this is a convergent sum, so E z is well defined. In fact, we show that E z is a field isomorphism. 
Note that this is the first use we make of the fact that the field of complex numbers C has characteristic zero. 



2.1. THE OPERATOR TOPOLOGY 
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Example 2.1.3 (Elementary D-Invariant Operator) For any pair of vectors a and (3 the elementary 
D-invariant operator from a to [3 — denoted Ep a — is the linear map on the logarithmic algebra defined for all 
a and 7 by 



where a is not a negative integer if a = (0). In the discrete case, the elementary D— invariant operators are 
continuous. 

Note that proj Q = E aa is the projection map T —* 2 a . In other words, 



These projections commute with D a . Note however that not all continuous, linear projections which commute 
with D a are expressible in terms of these projections. 

Continuous A subspace J of I + is said to be D-invariant if it is invariant under the fractional derivative 
D a for all real numbers a. An operator 9 on a D-invariant subspace J is said to be D-invariant when 
D a 9 = 9D a for all real a. 

Discrete A subspace J of T is said to be D-invariant or shift-invariant if it is invariant under the derivative 
D, or equivalently if it is invariant under E z for all complex numbers z. An operator 9 on a D-invariant 
subspace J is said to be D-invariant or shift-invariant when D9 = 9D, or equivalently if DE Z = E Z D 
for all complex numbers z. 

For example, all the operators mentioned above are D-invariant. 

A continuous linear operator 9 on T or I + is said to be a regular operator if it commutes with every 
elementary D-invariant operator E a p (except possibly when (3 = 0), that is, such that 



for [3 0. For example, D a and E z are regular operators.; 

A regular D-invariant operator on 2 + is called a Artinian operator, and one on all of T is called a 
differential operator. The set of Artinian operators is denoted by A + , and the set of differential operators 
by A. Beware that this notation is inconsistent with the notation for the corresponding concept in [24]. The 
notation used here was chosen because it is more logical. 

Clearly every differential operator restricts to an Artinian operator. Thus, A+ C A 

Proposition 2.1.2 1. The set of Artinian operators A+ is a complete topological ring in the operator 
topology ofI + . 





0E a p = E af3 9 
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2. The set of differential operators A is a complete topological ring in the operator topology of I. 

Proof: A+ and A clearly are AT-algebras, so it suffices to show that the limit of any Cauchy sequence of 
Artinian (resp. differential) operators is again an Artinian (resp. differential) operator. 

Let (9 n ) n >o be such a Cauchy sequence, and let 9 be its limit. Now, 

9E z p(x) = lim 9 n E z p{x) 

n—t+OG 

lim E z 9 n p(x) 

n—t+oo 

= E z lim 9 n p(x) 

n — >+oo 

since E z is a continuous operator. This in turn equals E z 9p(x), so 9 is D-invariant. Mutatis mutandis, we 
have E st 9 — 9E st , so 9 is regular. 1 

Our objective is to obtain structural characterizations of Artinian and differential operators. 



2.2 Taylor's Formula 

We shall derive analogs of Taylor's formula in the Logarithmic algebra 1. We begin by giving the following 
alternate definition of the shift operator: 



Proposition 2.2.1 For all complex numbers z, E z is a well defined continuous field isomorphism of I 
which fixes all constants. . 

Proof: Need only check that E z (p(x)q(x)) = {E z p(x))(E z q(x)) for all p(x), q(x) G I. 
E z (p(x) q (x)) = J2^V k (p(x) q (x)) 

k>0 

= E (*) (D»rt*)) (D-tf*)) 

fe>0 n+m=k ^ ' 
n+m 

= E (d™«(*)) 

n,m>0 

- (E z p(x))(E z q(x)). I 



2.2. TAYLOR'S FORMULA 
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Proposition 2.2.1 shows that E z satisfies the conditions of the characterization of Artinian 
composition in [25]. Any continuous field isomorphism of I which fixes all constants is determined 
by its values on the iterated logarithms Thus, E z is the O-composition associated with the 
the substitution of x + z for x, and logx + ^Zj >0 {—^y +1 z^ / jx J \ and so on. In general, 



z — ' n\ 



n>0 



= E 



(0,... ,0,1,0,...) 



n>0 



[-n]\n\ ~ n 



(x) 



4-E E 



n>0 t(p)=k 
Pk = l 



-z) n s{-n,pi) 
\n\ In 



k-2 



Y[e Pj - Pj+1 (-l,...,- Pj + l) 

3 = 1 



k-2 
3=1 



'Pi ~Pi 



= 4+ E (-i) pl+|p| 

i(p)=k 

Pk=l 

= 4+ E (-i)" fc - 1+|p| 
l(p)=k 

Pk=l 

where p* = (pi, ... , Pk-i)- End of Digression. 



+1 (-l,...,- Pj + l) 



k-2 

Yle^-p^il^.^pj-l) 

3 = 1 



g(-n),p* 

z n (n- l)\s(-n, Pl ) g( _ n) ^. 
o " 

z n (n- l)!s(-n,pi) , ( _ n) 



n>0 



E 

n>0 



:^(-™)>p* 



n 



Observe that E Zl E Z2 = e ZlD e Z2D = e {zi+Z2 ^ = E Zl+Z2 . 



Discrete By applying the shift operator to the harmonic logarithms of order (1) and degree n a nonnegative 
integer, \„\x), we obtain the following identity: 



(1 +„)" (Iog(l + a)-l-l l - 



[x + a) n [ log(x + a) - 1 - - - 



n 



a n - l x l [\ogx- 1- - 



i>n 
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and therefore: 

(l + o) n Iog(l + o) 



((1 + ar-l) 



L-n]! 



E 



L-ili 



+ 



EL 



= ((1 + a)"-!) 1 + - + ... + - -na 1 + - + ••• + 



n 



1 



7i 

2, 

n+l 



a' 1 + 



1 



n 


+ a"+ 2 


n 


-1 




-2 



n- 2 
+ • 



2 

3 /n 
2 \ 2 



n- 1 



a"- 2 - na"" 1 



Contrast with [35, Corollary 5.6]. 



2.3 Real Analysis 

Proposition 2.2.1 can be summarized as stating that 

E z p(x) = e zB p(x) = p(x + z) 

for all logarithmic series p(x) G X. Note that this identity is tautological, since we cannot "evaluate" the 
variable x. However, in the case of complex numbers, we can, and we have 

Proposition 2.3.1 All formal power series of logarithmic typep(x) € T represent asymptotic expansions 
of a complex valued function p(x) in a neighborhood of infinity. Moreover, for all real numbers c, E c p(x) 
represents an asymptotic expansion of p(x + c). 

Conversely, any asymptotic expansion of a function in a neighborhood of infinity relative to the ladder of 
comparison given by X"(x) is a formal power series of logarithmic type. 



Proof: The asymptotic expansions follow from Proposition ??. We conclude by observing 

E c p(x)\ x=z = e cD p(x) 



fe>0 



fe>0 



2.4. CHARACTERIZATION OF VARIOUS CLASSES OF OPERATORS 

2.4 Characterization of Various Classes of Operators 
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2.4.1 Artinian Operators 

Every Artinian operator 9 maps X a into itself for every nonzero vector a. We denote by 9 its restriction to 
I a , by an abuse of notation, and we say that 9 is an Artinian operator of X a into itself. 

As planned, we can now characterize the algebra of Artinian operators and its topology. It is isomorphic 
to an algebra (mentioned in [25]) we denote C < x >; it is in a sense the dual of the Noetherian algebra of 
Definition ??. 

Discrete In the discrete case, K < x > is the algebra of Laurent series in the variable x with complex 
coefficients. 

Continuous In the continuous case, K < x > is the algebra of Artinian series in the variable x with complex 
coefficients. 

_ jJ2 aeR c a x a : c a £ K and for all a £ R, there exists finitely manyl 
<2;> ~\&>0with Ca ^ Q j- 

Its operations and topology are similar to that of the Noetherian algebra. 

Theorem 2.4.1 The algebra of Artinian operators A + is naturally isomorphic to the Artinian algebra 
C < x > as a topological algebra A+ = C(D) Ro . 

Proof: (C(D) Ro C A+) Note that D a D b = D b D a , and B a E a/} = E afj D a . 

(A+ C C(D) R °) Let 9 £ A + be an Artinian operator. By regularity, 9 is determined by its action on the 
harmonic logarithms of order a for any particular (0). 

6\ a a {x)=Y,ea b \?{x). 
b 

Notice, that for any a and 6, e a b' only for finitely many b' > b. Next, 

hl /.,.A;'.:.n 



[b-a]\ 

D a 9\?(x) = D a ]Te bc A«(x) 
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Thus, 



e a -b, c = [ C ^ a ] e b,c+b In particular, setting a = b, 



eo,c 



e-b,c+b- 

Hence, 9 is determined by the eo, c , and therefore equals the Artinian series J2 C e c,oD c / L c l'- 

(Continuity) It remains only to show that Artinian series in the derivative are continuous in the operator 
topology (§2.1), and that Artinian operators are continuous in the Artinian topology. However, c„D a ™ 
converges whenever a n and c„ converge, and c„D a ™ converges to zero whenever a n increases without bound. 
Moreover, all Cauchy sequences of operators are linear combinations of the above. I 

If /(D) is an Artinian operator which is represented by a delta series (that is, a series of degree one), 
then /(D) is called a delta operator. 

Discrete Recall the important fact that if /(D) is a delta operator, then the sequence of powers (/(D) n )„ e z 
is a pseudobasis for the field of Artinian operators A + . That is, for every Artinian operator g(D) € A + 
of degree d there is a sequence (ck)k<d indexed by integers k < d such that 5(D) = J2k<d c fc/(D) fc . 
This sequence can actually be calculated via Theorem 3.3.3. The notation /(D) a; ™ comes from the 
continuous case (see below); in the context of the discrete theory, it should be read simply as /(D) a . 

Continuous Recall the important fact from [25] that if /(D) = J2 a eR CaXa 1S an Artinian operator of 
degree d 7^ (for example, a delta operator), and n is an integer, then the set {/(D) o; ™ : a £ R} is a 
pseudobasis for the field of Artinian operators A + where 




\-d a 
— X 



M x " 7 \aeM Cd J 

and m = \cd\ and — arg(c<i) are the modulus and argument of the leading coefficient of /(D). 

That is, for every Artinian operator g(D) € A + there is a sequence (c a ) a(E R such that ^2 aeR c a/(D) a; " 
converges to g(D). When /(D) is a delta operator, this sequence can be calculated via Theorem 3.3.3 
or the methods of [20]. 

Similarly, the series of nonnegative integers powers of any delta operator /(D) is a pseudobasis for A. 



2.4.2 Differential Operators 

In analogy with the preceding result for Artinian operators, we obtain the following structure theorem for 
differential operators: 

COROLLARY 2.4.2 The topological algebra of differential operators is naturally isomorphic to the topo- 
logical algebra of formal power series in the derivative as a topological algebra: A = C[[D]]. 

Proof: Such series are the only members of A + which are well defined on T^ . !S 

As opposed to Artinian operators which are invertible if nonzero, a differential operator is invertible in 
A if and only if it is of degree 0. 



2.4. CHARACTERIZATION OF VARIOUS CLASSES OF OPERATORS 



15 



2.4.3A D-invariant Operators 



The rings of continuous, linear, D-invariant operators (which are not necessarily regular) over 1 and 1+ 
are structured as follows: 

Proposition 2.4.3 LetlZ andlZo be the rings of D-invariant linear operators onl + andl respectively 
which are continuous on each l a . Then: 

Discrete 1. 1Z is the closure D, in the operator topology, of the span of the operators D n E a p where a and 
(3 are nonzero vectors with finite support of integers, and n is an integer. 

2. TZq is the closure T>q, in the operator topology, of the span of the operators Y) n E a fj where 

(a) a and (3 are vectors with finite support of integers, 

(b) n is an integers, 

(c) a ^ (0) unless [3 — (0), and 

(d) [3 ^ (0) unless n is a nonnegative integer. 

Continuous 1. 1Z is the closure T>, in the operator topology, of the span of the operators T} a E a /3 where a 
and (3 are nonzero vectors with finite support of real numbers, and a is a real number. 

2. IZq is the closure T> , in the operator topology, of the span of the operators D a E a p 

(a) a and [3 are vectors with finite support of real numbers, 

(b) a is a real number, 

(c) a 7^ (0) unless (3 = (0), and 

(d) [3 ^ (0) unless a is a nonnegative integer. 

Proof: (DCR and T> C TZq) Elementary D-invariant operators commute with the derivative. Thus, 
they commute with all Artinian and differential operators. Hence, they are in fact D-invariant. 

Observe that E a pD a is continuous, linear and D-invariant for a, (3 ^ (0). We conclude that every operator 
in V and T> is continuous, linear, and D-invariant. 

(TZ C V) Let 8 be a continuous, linear, D-invariant operator. For each pair of vectors a, [3 ^ (0), define 
S a p = proj^^proj^. Obviously, 9 = ^(o) @af3- It suffices to show that for all nonzero vectors a and [3, 
there is an Artinian operator / a /3(D) e A + such that 9 a p = f(D)Ep a . 

However, Ep a 6 a p is a continuous, linear, D-invariant operator on 2 a , so Ep a 9 a/ 3 = /(D)proj ( g for some 
Artinian operator /(D). Hence, 6p a = Ep a f(T))E aa = f(D)Efj a as desired. 

(TZo C T> ) Similarly, it suffices to show that 9 a .(a) (as defined above) is equal to zero for a ^ (0). Assume 
not towards contradiction. By the reasoning above, Qi (o) = /(D)-E a .( ) for some nonzero differential operator 
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/(D). Let g(D) be the inverse of /(D). Since g(D) is a D-invariant operator, the product g(D)0 is also D- 
invariant. Hence, without loss of generality, # a> (o) = ^ a ,(o)- We calculate that _E a> ( )Dl — E a ^ )0. However, 
we also know that DE a ( )1 = Ei^' a ^ 0. Contradiction. I 

We omit a discussion of non-D-invariant continuous, linear, regular operators and their expansions in 
terms of D and a (to be defined later) , since this would be an unnecessary digression. 

However, there are other interested related questions left unsolved. 

Open Problem 2.4.4 Is there a simple characterization of non-linear, continuous, shift-invariant op- 
erators on polynomials? Or of the logarithmic algebra? For example, p(x) ^ p{x) 2 . 



2.4.4A Differential Equations 



We next make some general remarks about solutions of differential equations of infinite order (and thus, 
in particular, of difference equations). We have seen (Theorem 2.4.1) that the algebra of Artinian operators 
is isomorphic to the Artinian algebra in the "variable" D as a topological algebra over the complex numbers. 
In particular, every Artinian operator is invertible. Hence, every differential equation of the form 

f(B)p(x) = q(x) (2.1) 

where /(D) G A + is an Artinian operator and q(x) € X + has a unique solution p{x) in X + . For example, 
q(x) may be any rational function of x whose numerator is of smaller degree than the denominator. 

If q(x) e X (for example, if q(x) is an arbitrary rational function) the solution may not be unique. We 
can nevertheless define the natural solution of equation (2.1) as follows. Let n be any bijection between 
the equipotent sets A and A — {(0)} where A is the set of all vectors a with finite support of reals (resp. 
integers). Define P v to be the linear map (continuous on each X a ) defined by 

P-rr = ^ ^ E-xa,a 
a 

or equivalently 

PMx) = \™(x) 

For example, we might have 

Pit = X] ^(n)(n+l) + X/ Eaa 
n>0 ajt(n) 

where the second sum in each equation is over vectors which do not consist of a single nonnegative integer. 
Thus, Ptt : X — > X + and its inverse P^T 1 : X + — > X is given by P~ 1 X"(x) = a {x). Now, let s(x) = P 7T q(x), 
and consider the differential equation 

f(D)r(x) = s(x). 



2.5. AUGMENTATION 
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Since s(x) G I + , this differential equation has a unique solution r(x) G J + . Now, set p(a;) = P~ 1 r(x) 
to obtain a solution of equation (2.1). The present definition agrees with (and is simpler than) all other 
definitions given of a natural solution over the complex numbers. Moreover, since p(x) does not depend on 
the choice of n, the solution has been chosen naturally. 

A notable example is the difference equation 

Ap(x) = l/x. 

By the above remarks, it has a unique solution in X + , which turns out to be the V'-function ip(x), the 
logarithmic derivative of the gamma function. Thus, the theory of the ^-function can be developed purely 
formally. (See chapter 5.1.1.) 

So far, we have only consider linear differential equations. This leads us to the following open problem. 

Open Problem 2.4.5 In general, what differential equations have solutions? And when do they have 
canonical solutions? 

The logarithmic algebra is not a differentially closed field since Dp(x) = p(x) only has one solution — 
p(x) = 0. However, we can redefine the degree of a differential operator so that we are consistent with [25] by 
insisting that its degree is the lowest (rather than highest) exponent of D with a nonzero coefficient. Under 
these circumstances, Dp(x) — p(x) is no longer a counterexample. 

Open Problem 2.4.6C Under this definition of degree, is the logarithmic algebra differentially closed? 
And if so, can we modify the logarithmic algebra so that it is differentially closed under the usual definition? 

Recall that differentially closed fields are known to exist in every character; however, no example of a 
differentially closed field of character zero has been found. 

2.5 Augmentation 



Although it is impossible to evaluate at zero any expression involving logarithms, the following definition 
nevertheless serves as the logarithmic analog of evaluation at zero. 

Definition 2.5.1 (Augmentation) For each vector a, we define the augmentation of order a to be the 
linear functional (} a (continuous on eachX a ) from the logarithmic algebra to the complex numbers such that 

or equivalently using the notation introduced after Corollary ??, 



(p(x)) a = [^(x)]p(x). 
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We digress to indicate the algebraic significance of augmentation in the discrete case. First, 
note that in this case () a is continuous. 

When () is restricted to , the augmentation reduces the evaluation of a polynomial at 
x = 0. That is, (p(x)) Q = p(0) for p(x) G 2^ . An augmentation of nonzero order can be viewed 
as a generalization of evaluation at x — 0; it is closely related to the residue of complex variable 
theory. (Recall that the residue of a Laurent series is its coefficient of a; -1 .) For instance, for 

(p(x)\ = Res(Dp(aO). 
In fact, for a a vector of j nonnegative integers and p(x) e l a , 




Note that equation (2.2) also holds for all p(x) e /3<ct 2^ /3 . End of Digression. 



We derive formulas relating the augmentation to the derivative. 



Proposition 2.5.2 1. For all a and b, and all a, (3 ^ (0) 

(V b \ a a (x)) fj = [a\\8 ab 8 aP . 

2. For all a and all nonnegative integers n, 

3. If /(D) is an Artinian operator with coefficients [D a ]/(D) = c a , and p(x) is a formal power series of 
logarithmic type with coefficients 

b a a = [\™(x)]p(x), 
then the augmentation (f(D)p(x)) a is given by the finite sum 

(f(D)p(x)) a =J2laV-c a bZ. 

a 

4- ///(D) is a differential operator with coefficients Ck = [D fc ]/(D), and p{x) is a formal power series of 
logarithmic type with coefficients 

b: = [\:(x)] P (x), 

then the augmentation (f(D)p(x)) a is given by the finite sum 

(f(D)p(x)) a =J2klc k b a k . ■ 

fc>0 
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Some special cases are of interest. The augmentation of the derivative of a formal power series of 
logarithmic type can be described by (D a p(x)) a — [a]lb" where p(x) is as in Proposition 2.5.2. Similarly 
the augmentation of the action of a differential operator on a harmonic logarithm is (/(D)A"(x)) a = La~|!c„ 
where /(D) is as in Proposition 2.5.2. 

The augmentation leads us to a version of Taylor's formula for formal power series of logarithmic type: 
Theorem 2.5.3 (Taylor's Formula) Let p(x) € 2 + . Then we have the following convergent expansion in 

Proof: By linearity and continuity it suffices to consider the case p(x) = A" (a;). This special case was 
treated above (Proposition 2.5.2). ffl 

Note that the series of the right hand side is convergent. 



Discrete Observe that in the discrete case, any formal power series of logarithmic type is determined by 
the augmentations of its derivative. 

Continuous Nevertheless, in the continuous case, members of can not be recovered from the augmen- 
tations of their derivatives. For example, (f(D)y/x)^ = for all differential operators /(D). 

On the other hand, in cither case, an Artinian operator can be recovered from the augmentations of 
its action on the harmonic logarithms of any particular order a ^ (0) and a differential operator can be 
recovered from the augmentations of the harmonic logarithms of any particular order as is demonstrated by 
the following theorem. 

Theorem 2.5.4 (Expansion Theorem) 

1. Let /(D) be an Artinian operator, and let a ^ (0), then we have the following convergent expansion 



f(D) = V Ui J "^ p" 
n ' ^ \a]\ 



2. Similarly, if /(D) is a differential operator, and a is a vector, then we have the following convergent 

S&TZ6S 

/(D) = E iZMSlk D n , 



n>0 
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The following argument is used repeatedly in the next two sections — often implicitly. 
PROPOSITION 2.5.5 (Spanning Argument) 

1. Let p(x) e I + . If (f(D)p(x)) a = for all vectors a ^ (0) and all Artinian operators /(D), then 
p{x) = 0. 

2. Let a ^ (0) be a vector, and let /(D) be an Artinian operator. If (f(D)p(x)) a = for all p(x) £ l a , 
then /(D) =0. 

3. Discrete Similarly, in the discrete case, for let p(x) G l^K If (f(D)p(x)}^ = 0, for all 

differential operators /(D) e A, then p(x) = 0. 

4- Let a be a vector, and let /(D) be a differential operator. If for all p(x) G I a , (f(D)p(x)) a = 0, then 
/(D) -0. I 

Open Problem 2.5.6 Is there a simple formula expressing any monomial i a or product of harmonic 
logarithms X"(x)X^(x) in terms of harmonic logarithms? That is, is there a simple way to calculate the 
coefficients of l a or \"(x)\%(x) given by Theorem 2.5.3? 

Discrete Note that the expansion of \^n \x)Xm\x) is known. 



Chapter 3 



Roman Graded Sequences 



3.1 Graded Sequences 



This chapter is devoted to the study of the logarithmic analog of sequences of polynomials of binomial 
type. However, before one can walk; one must crawl. We must first define the logarithmic analog of a 
sequence of polynomials. Classically, a sequence of polynomials (p n {x)) n >o is subject to the requirements 
that deg(p n (x)) = n for all n. Here the requirements are slightly more complicated, yet the idea remains the 
same. 

Definition 3.1.1 (Graded Sequences of Logarithmic Series) The sequence : a real (resp. 

an integer) and a a vector of reals (resp. integers)} is called a continuous (resp. discrete) graded sequence 
of formal power series of logarithmic type if the following six (resp. five) conditions hold. (Part of a typical 
such sequence is illustrated by Figure 3.1 

1. For all a and a / (0), p"(x) is a homogeneous formal power series of logarithmic type of order a. For 
example, the series in the k^ 1 row of Figure 3.1 are homogeneous of order (k). 

2. For n a negative integer, Pn\x) = 0. That is, the left side of the middle row of Figure 3.1 is filled 
with zeroes. 
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Figure 3.1: Part of a Typical Graded Sequence of Formal Power Series of Logarithmic Type 



... pL" 3 2) (z) 




p^Hx) 


p£" 2) (x) 


p5~ 2) w 


p(" 2) (x) 


p ( 3 - 2) (x) ■■■ 











Po(x) 


Pl(x) 


p 2 (a;) 


p 3 (ir) 


p-3(x) 


p-2(x) 


p-i(x) 


p^x) 


p[ 1} (x) 




pW(x) ••■ 


■■ P ( *l(x) 


pL 2 >(x) 


P L a i(x) 


p< 2 >(x) 


p< 2 >(x) 


pf\x) 


pf(x) ••■ 



5. For a not a negative integer, p a °\x) is a logarithmic series of degree a. In other words, in the discrete 
case for n a nonnegative integer Pn (x) is a polynomial of degree n. That is, the right side of the 
middle row of Figure 3.1 is filled with polynomials in order of degree. 

4. For a 7^ (0), p"{ x ) is of degree a. For example, the n^ 1 column of Figure 3.1 consists of series of 
degree n. 

5. (Regularity) For all a and 13, and all a 7^ (0), E a ppP(x) —p"{x). For example, the operator i%)(fc) 
sends the k^ 1 row of Figure 3.1 to the row. 

6. Continuous (Continuity) In the continuous case only, the map from the reals to I a defined by 

!p"{x) a not a negative integer 
(3.1) 
Pa a is a negative integer 

must be continuous for all a. 

When a sequence p" (x) has the property that its leading coefficients are all positive real numbers, then it is 
lid to be standard. 



The logarithmic series p_\(x) is called the residual series of the graded sequence p^{x). (Indicated by a 
box in Figure 3.1.) 



The principle subsequence of p"(x) is the subsequence {Pa(x)) n£l{ defined by 

pi°^ (x) for a not a negative integer, and 
pi^ (x) for a a negative integer. 



Pa(x) 



For example, 
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Proposition 3.1.2 A"(x) is a standard graded sequence of formal power series of logarithmic type. Its 
residual series is 1/x, and its principal subsequence is the sequence of powers of x, {x n ) neZ - 

Proof: The only thing worth checking is that the sequence is continuous. However, since the Gamma 
function is infinitely differentiable at all points other than the nonnegative integers, the map a \— ► s(a, k) is 
continuous. Hence, the map defined by equation (3.1) is continuous. I 

The nonzero elements of any graded sequence form a pscudobasis for I, and the restriction of a graded 
sequence to any particular a forms a pseudobasis for X a . Thus, for any pair of graded sequences p" (x) and 
q" (x), there is a unique continuous linear operator 9 such that Op"(x) — q% (x). We study such operators in 
detail in chapter 4. 

Theorem 3.1.3 Letp"(x) be a graded sequence. Then every logarithmic series g(x) e I can be uniquely 
written as an expression 

g{x) =J]5XW 

a, a 

where g a (x) G T a . This is a convergent expansion in the topology of Noetherian series, and an asymptotic 
expansion in the topology of the complex numbers as x tends towards infinity. ffi 

We are able to actually calculate the constants b™ in many cases using such results as Theorems 3.3.4 
and the corresponding result in [22]. 



3.2 Roman Graded Sequences 

In this section, we introduce the central concept of this work. It is known that the operational calculus of 
formal differential operators is intimately associated with sequences of polynomials of binomial type, that 
is, with sequences of polynomials p n (x) satisfying the binomial identity (Definition 3.6.1) 

p n (x + a) = ^2 ( 2jPk(x)p n -k(a) 

A good many sequences of polynomials occurring in combinatorics and in the theory of special functions turn 
out to be of binomial type. For example, the powers x n , the lower factorial (x)„, the upper factorial (x) n , 
the Abel polynomials A n (x), the LaGuerre polynomials L n (x), and the inverse- Abel polynomials fj, n (x) are 
all sequences of polynomials of binomial type. We give here the logarithmic generalization of this notion; 
such graded sequences of formal power series of logarithmic type are called Roman graded sequences. We 
derive five equivalent characterizations of such graded sequences. We anticipate the fact (Theorem 3.6.4) 
that the five notions introduced below coincide: 

1. Roman graded sequence (Definition 3.2.2), 
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2. associated graded sequence (Definition 3.3.2), 

3. basic graded sequence (Definition 3.4.1), and 

4. conjugate graded sequence (Definition 3.5.1), 

5. graded sequence of logarithmic binomial type. (Definition 3.6.2) 

Wee motivate the definition of a Roman graded sequence, by deriving a formula for the action of a product 
of two Artinian operators on the harmonic logarithm. 

Proposition 3.2.1 1. Let /(D) and g(D) be Artinian operators. Then for all a and for all a ^ (0) 
we have the following finite sum: 

(f(D)g(D)>£(x)) a = ]T [I] </(D)A?(aO> Q (g(D)\" a _ b (x)) a . (3.2) 

b 

2. Similarly, let /(D) and g(D) be differential operators. Then equation (3.2) holds for all vectors a. 
Proof: Let c a = [D a ]/(D), and d a = [D a ]g(D). The product of the two series is given by the sum 



f(D)g(D)=J2^Ecada-b^ B b . 



Hence, by Theorem 2.5.2, 

(f{D) 9 {D)\?{x)) a = [bV.^Cadb-a 



b 



E 



(f(I>)>£{x)) a (g{D)\?_ a (x)) a 



The extension to products of more than two operators follows easily by induction. 

We introduce Roman graded sequences by the following definition. It will shortly be seen that simpler 
alternate definitions can be given. 

Definition 3.2.2 (Roman Graded Sequences) Letp"{x) be a graded sequence of formal power series of 
logarithmic type. The graded sequence is a Roman graded sequence if for all a and a, we have the following 
finite sum 

(/(D) ff (D)p°(x)) a = £ [I] (f(D)Pb(x))a (9(VK_ h (x)) a (3-3) 

b 

where /(D) and g(D) are Artinian operators if a ^ (0) ; and are differential operators if a — (0). 
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For example, A" (a;) is a Roman graded sequence. 

Proposition 3.2.3 A graded sequence is Roman if and only if equation (3.3) holds when /(D) = D a 
and g(D) = D b . 

Proof: Linearity and continuity. I 

Proposition 3.2.4 A graded sequence p"(x) with coefficients d a b — [A" {x)]p" (x) is Roman if and only 
if for all a, b, and c: 

d c ,a+b = ^2 debd c ~e,a. (3.4) 

e 

Proof: We demonstrate the "only if;" the reasoning for the other implication is similar. 

(D a D b p«(x)) a = [a + b]\d c , a+b 

= [a-\\[b-]\J2[ C e ]<ad c -e,b 
e 

= E^l< D ^)U D ^-eM) Q - 

e 

Note that in the case a = (0), we need only consider nonnegative integers a and b, so the same argument 
applies. 1 



3.3 Associated Graded Sequences 



We proceed to derive an altogether different characteristic property of Roman graded sequences, which 
uses delta operators. Such differential operators may be viewed as playing the role of the derivative — much 
like the forward difference operator A (Definition 5.1.1) acts on the polynomial sequence of lower factorials 
(x) n = x(x - 1) • • • (x - n + 1). 

Proposition 3.3.1 Let /(D) (and n an integer) be a delta operator in A + . Then there is a unique 
graded sequence of formal power series of logarithmic type, p"{x), such that for all a, b and a 



a + b 
a 



( Discrete (f(D) a p%(x)) a = la]\5 ab | 
\ Continuous (/(D) a ->«(x)) Q - [a]\5 ab . J 



(3.5) 
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Proof: (Uniqueness) Spanning argument (Proposition 2.5.5.) 

(Existence) We need only consider the continuous case. Let 4 a) = [D b ]/(D) a;rl be the coefficients of 
/ (D) a;n , and let c b = = [D b ]f (D). Thus, c{ a) ^ for all real numbers a. Similarly, let d% b = [\%(x)]p%(x) 
be the coefficients of p" (x). In this notation, equation (3.5) is equivalent to the equation 

b 

Hence, the coefficients d" b can be computed recursively as: 

d "» = TTTTTbT \^» - E ^^dab J • 

l b V-al \ ob J 

The recursion is well defined since d" h is a Noetherian sequence in b. Finally, notice that (x) is regular by 
symmetry, and that deg(p^(x)) — n since b" n = [n] \/\n\ lail 1 ^ ^ 0. Finally, note that the map defined by 
equation (3.1) is continuous since a ^ /(D) a; ™ is continuous. 1 



Definition 3.3.2 (Associated Graded Sequence) Let /(D) be a delta operator (and n an integer). The 
unique graded sequence mentioned in Proposition 3.3.1 is called the (n^h) associated graded sequence of the 
delta operator /(D). We also say that the sequence is n— associated with /(D). (The term ri^ 1 applies only 
in the continuous case.) 



For example, by Proposition 2.5.2, \"{x) is the standard graded sequence associated with the delta 
operator D. 

We can now generalize Theorem 2.5.4 to explicitly determine the coefficients in the expansion of an 
arbitrary Artinian operator in terms of the powers of a delta operator: 



Theorem 3.3.3 (Expansion Theorem) Let the graded sequence p"{x) be n-associated with the delta 
operator /(D). Then for all Artinian operators g(D) and vectors a ^ (0) we have the following convergent 
sum. 



Discrete g(D) 
Continuous g(D) 



-f(B) k 



J2 (9(D)Pa(x)) a j^a:n 



(3.6) 



When a = (0), equation (3.6) holds for all differential operators <?(D). 
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Proof: By Definition 2.5.1 we have, 



E 



La]! 



/(D)°'»p?W 



- E 



(sW(4 



(/(D)^"K(a;)) Q 



La]! 

= (ff(D)p°(x)) a . 

The conclusion follows by the spanning argument. IH 

Dually, we obtain the explicit form of the expansion of an arbitrary formal series of logarithmic type 
as a linear combination of elements of a Roman graded sequence. This gives a useful generalization of 
Theorem 2.5.3. 

Theorem 3.3.4 (Logarithmic Taylor's Theorem) Letp"(x) be the (rfi 1 ) graded sequence associated with 
the delta operator /(D). Then for every formal power series of logarithmic type p(x) e 1 + we have the 
following convergent sum 

Discrete p(x) = ^ E^^^) 



Continuous p(x) 



ct^(O) k 

EE 



[k]\ 



Proof: We need only consider the continuous case, but first a simple lemma. For all p(x), q(x) el + , by 
Part (1) of the spanning argument (Proposition 2.5.5), it is clear that p{x) = q(x) if and only if for all a 
anda^(O), (D a p(x)) a = (D a q(x)) a . 



When we apply the Expansion Theorem where 5(D) is set equal to D a , we find that 

= E (Da ff L /( D ) b; " 



Thus, 



(D>(x)) Q 



E 

b 



L6]! 



- (f(D) b -- n p(x)) a 
(/(D)^>(x)) a ' 



Therefore, the above remarks, we conclude that 



In the continuous case, where this is relevant, we may now classify standard associated sequences (Defi- 
nition 3.1.1). 
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Proposition 3.3.5C The nft 1 associated sequence of a delta operator /(D) is standard if and only if the 
leading coefficient of /(D) is a positive real number and n = 0. § 

Such delta operators are be called standard operators. 



3.4 Basic Graded Sequences 



Definition 3.4.1 (Basic Graded Sequence) Let /(D) be a delta operator. A graded sequence of formal 
power series of logarithmic type p"(x) is called the (n tn ) basic graded sequence for /(D) if for all a, 

1. (p%(x)) a = 1: 

2. For all a > (and thus all a ^ 0), (p™(a;)) Q = 0; and 

3. Discrete For all integers n, and all vectors a of integers with finite support, f(D)p^(x) — \n\ p"_ 1 {x). 

Continuous For all real numbers a and b, and all nonzero vectors a of real numbers with finite 
support, 

f(D) b -<"p« a (x) = J^ P Z_ b { x ). (3.7) 

Theorem 3.4.2 1. Let p"(x) be a logarithmic graded sequence. Such a graded sequence is the (rfi 1 ) 
basic graded sequence for the delta operator /(D) if and only if it is the (vft 1 ) associated graded sequence 
for the delta operator /(D). 

2. Every delta operator has a unique (rfi 1 ) basic sequence. 

3. Every basic sequence is basic for a unique delta operator (and integer n). 

Proof: (1: If) Properties 1 and 2 of Definition 3.4.1 follow from Definition 3.3.2. Property 3 follows 
from the following series of equalities, and the spanning argument as shown here in the continuous case: 

(f(D) ain {f(V) b '' n Pc(x)))a = </(D) a+b; >«(x)) Q 

= lc\l5 c ,a+b 

= T^]T </( D ) a; >c-6(*)> a 
(1: Only if) We proceed differently in the continuous and discrete cases: 
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Discrete By induction we have, f(D) k p"(x) — [n]\p"_ k (x)/[n — k~\l. Hence, 

</(d)*k(z)> q = (tS^pn-^y 

= [n\\5 n ^ k . 

Hence, p„(x) is the associated graded sequence of /(D) by definition. 
Continuous By equation (3.7), (f(Df' n p a a {x)) a = j^S a _ bfi = L«l^a6- 

(2 and 3) Immediate from 1. S 

The simplest example of a basic graded sequence is the graded sequence of harmonic logarithms A" (a;); 
by Theorem ??, it is the standard basic graded sequence for the delta operator D. It can be viewed as the 
natural logarithmic extension of the sequence of powers of x. More generally, every sequence of binomial 
type (and every factor sequence) has a natural logarithmic extension into a basic graded sequence related to 
the same delta operator. 

D is invertible on I + . Hence, just as in the continuous case (equation (3.7)), we have: 



PORISM 3. 4. 3D Let p"(x) be the associated graded sequence of formal power series of logarithmic type 
of the delta operator /(D), then for all integers n and k, and all nonzero vectors a with finite support of 
integers. 

f(D) k p«(x) = J^p^W (3.8) 
Moreover, when k is a nonnegative integer, equation (3.8) holds for alia. B 



We may now connect the notion of an associated graded sequence with that of a Roman graded sequence. 



Proposition 3.4.4 Let p"(x) be an associated graded sequence, and let /(D) e A+ be an Artinian 
operator. Then we have the following convergent sum 

f(D)p a a (x) = J2 [l] </(DK(*M-6(z) 

b 

for all a ^ (0) ; and also for a = (0) if /(D) is in fact a differential operator. 
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Proof: We need only consider the continuous case, suppose p% (x) is n-associated with the delta operator 
g(D). Then for all b 

= E [cl (^ D ) 6; >cW) a K-cW- 

c 

Since {g(D) 6;n : 6 £ R} is a pseudobasis, we can use continuity and linearity to replace g(D) b;n by 
/(D). I 

The following identity is classically proven [40] by introducing a new variable. However, we have a much 
simpler proof. 

COROLLARY 3.4.5 Letp^(x) be an associated graded sequence with coefficients 

c ab = [XUxM(x). 

Then for all a, b, and a we have 

D^(z) = E^K- d (*). 

Proof: Proposition 3.4.4. li 

As another consequence of Proposition 3.4.4 we obtain: 

Proposition 3.4.6 A graded sequence of formal power series of logarithmic type p"(x) is an associated 
graded sequence if and only if it is a Roman graded sequence. 

Proof: (Only if) Let /(D) and 5(D) be Artinian operators. By Proposition 3.4.4, 

/(DKW = E [I] </(DK(z)><X- 6 (*)- 

b 

Thus, 

</(D) 3 (DK(x)) Q = ^(D) ^ [ a b ] (f(D)p%(x)) a p%_ b (x) 

= EiIl(/( D KWUs( D )^w)o- 
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(If) Conversely, let p"{x) be a Roman graded sequence. We define a sequence of Artinian operators 
/b(D) by the relation 

(^(D)^))^^]!^. 

It suffices to show that 

Discrete / 6 (D) = /(D) 6 '™ 1 
Continuous / 6 (D) = /(D) b J 
for some delta operator /(D) (and some integer n). 

By the spanning argument, /&(D) is well defined. Now, 

(/ b (D)AW(x)) (1) =0 

for a < b, and 

(/,(D)Al 1) (x)) (1) jLO, 

hence deg(/j,(D)) = b. In particular, /i(D) is a delta operator. Since p"(x) is a Roman graded sequence, we 
infer that 

(/ b (D)/ c (D)rf)(x)) (1) = E^l<^( D )^ 1) (-))(i)</ C ( D >i 1 - d (-)>(i) 

d 

= [a]\5 atb+c 

By the spanning argument, /j(D)/j(D) = / i+J (D). Hence, 
Discrete By induction, / n (D) = fi(D) n . 

Continuous By the characterization of exponentiation in [25], there is an integer n such that /&(D) = 
/i(D) b '" for all n. ■ 

Thus, every Roman graded sequence is associated with a unique delta operator (and integer) and visa 
versa. 



3.5 Conjugate Graded Sequences 

Each delta operator (and integer) has another graded sequence associated with it: its conjugate graded 
sequence. We see that these sequences are Roman; however, they are not associated with the same delta 
operator for which they are conjugate. 
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Definition 3.5.1 (Conjugate Graded Sequence) Let /(D) e A+ be a delta operator. Its (rfi 1 ) conjugate 
graded sequence g"(x) is defined as 



Continuous g a (x) = ^ ^ — ^A^(x) 

b<a L I- 



Discrete q n (x) = ^ " T^n —hi*) 




for all a and a. 

Indeed, for all delta operators /(D) (and integers n), the graded sequence (x) as defined by equation 
(3.9) automatically meets the conditions of Definition 3.1.1. 

The canonical example of a conjugate graded sequence is the graded sequence of harmonic logarithms; it 
is the standard conjugate graded sequence of the delta operator D, since by Theorem 2.5.3, 



In the above example, the n tn conjugate graded sequence of a delta operator is a Roman graded sequence. 
This fact is true in general: 

Proposition 3.5.2 A graded sequence of formal power series of logarithmic type q%(x) is Roman if 
and only if it is the (vft 1 ) conjugate graded sequence of a delta operator. Moreover, the delta operator (and 
integer) are unique. 

Proof: (If) We need only consider the continuous case. Let <?"( x ) be the n^h conjugate graded sequence 
of the delta operator /(D). Now, by equation (3.9) the coefficients c a b of p"(x) are given by 



\ a a (x) = Y,(v b K(x)) a K(x)/lbV- 



b 



Cab = [xUx)K(x) 




[b]\ 



It suffices to show that the c a b satisfy equation (3.4). For any b\, 62 




(/(D) 




since 



A" (2) is a Roman graded sequence. The last expression equals 



bi+b 2 
bi 



c aM+b2- Hence, equation (3.4) 



holds. 
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(Only if) Conversely, suppose that q"(x) is a Roman graded sequence. Let c ab denote the coefficients 
of q%(x) 

c ab = [\?(xM(x). 

Define the Artinian operator /d(D) by 

(f d (D)\«(x)) a = [d]\c ad . 

By the spanning argument, this condition defines /d(D). As in the proof of Proposition 3.4.6 which employs 
a similar technique, it suffices to observe that /d(D) = f(D) d (rcsp. /(D) d;n ) for some delta operator /(D) 
(and integer n). 

In fact, we have: 

Thus, /i(D) is a delta operator. Now, 

(fb 1+b2 (Vmx)) a = [bi + b 2 V.b aM+b2 
by definition, and that expression equals 

d 

by equation (3.4). This in turn equals 

(/ 6l (D)/ 62 (D)A«(o ; )) Q 
since A" (x) is a Roman graded sequence. In other words, 

(f bl+b2 (D)\Z(x)) a = (f bl (D)f b2 (D)\«(x)) a , 

so by the spanning argument f i+ j(D) — /j(D)/j(D). 

Discrete Hence, by induction, /„(D) = .^(D)". 

Continuous Notice first that d i— > /d(D) is a continuous map. Hence, by the characterization of exponen- 
tiation in [25], there is an integer n such that /t(D) = /i(D) 6; ™ for all n. 

In other words, q"(x) is the (n^h) conjugate graded sequence of the delta operator /i(D) (and no 
other). I 

In the continuous case where this is relevant, we may now classify standard conjugate sequences. 

PROPOSITION 3.5.3C The vft 1 conjugate sequence of a delta operator ./(D) is standard if and only if 
/(D) is a standard operator and n = 0. 1 
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3.6 Graded Sequences of Logarithmic Binomial Type 



We motivate our final reformulation of the Roman condition by briefly recalling some facts about se- 
quences of polynomials of binomial type. 

Definition 3.6.1 (Polynomial Sequence of Binomial Type) A sequence (p n (x)) n>0 of polynomials is 
of binomial type if each polynomial p n (x) is of degree n, and for all field elements a, and all nonnegative 
integers n, 



It is a basic result of Umbral calculus (see for example [40]) that every sequence of binomial type is 
associated with a unique delta operator and vice versa. (p n (x)) n >o is the associated sequence of the delta 
operator /(D) if and only if f(D)p n (x) = np n -\{x) for positive n, and p n (0) — <5„,o for nonnegative n. 

In the present theory, the logarithmic analogs of sequences of binomial type are the graded sequences of 
binomial type defined below. As before, this definition turns out to be equivalent to that of a Roman graded 
sequence. 

Definition 3.6.2 (Graded Sequence of Logarithmic Binomial Type) A graded sequence of formal power 
series of logarithmic type p"(x) is of logarithmic binomial type if for all a and a, and all complex numbers 

z, J2b Lbl \^' Z Pb ) ' > ( x )/(n)P'a-b( x ) * s a convergent sum which equals E z p"(x). 




Discrete Note that in the discrete case (E z p£>{x)J (n) is merely p\ '{z). Thus, for a = (0), we obtain 

simply the definition of a sequence of polynomials of binomial type. For a = (1) and a a negative 
integer, we obtain factor sequences. [39] 

Thus, the notion of a graded sequence of logarithmic binomial type subsumes both the notion of a 
sequence of polynomials of binomial type, and the notion of a factor sequence. In view of the present 
theory, these older notions can be seen as obsolete. 

We next prove that, as promised, every graded sequence of logarithmic binomial type is a Roman graded 
sequence, and conversely. 

Theorem 3.6.3 A logarithmic graded sequence p"(x) is the (rfi 1 ) basic graded sequence for some delta 
operator /(D) (and integer n) if and only if it is a graded sequence of logarithmic binomial type. 
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Proof: (Only if) Proposition 3.4.4. 

(If) To prove that a graded sequence is basic, we must demonstrate each of the three properties enumer- 
ated in Definition 3.4.1. 

(1) By Definition 3.6.2, 
and 1 is a field, so we have 

(2) Whereas for a positive, we have 

o =p^( X ) -p^\ X ) - e l;i (** o) (4o)A(*). 

and the p^\x) form a pscudobasis for X^°\ so (p^\x)j^ = for b > 0. Thus, by regularity, 
for all a and a. 

(3) Define a sequence of continuous, linear operators Q b by the identity 

Q b Pa( x ) = J^PZ-bW- 

Clearly, Q C Q C = Q b+C , and the map b Q b is continuous in the operator topology, so by the characterization 
of exponentiation in [25], it remains only to show that Q — Q 1 is a delta operator. 

By inspection Q is regular, and lowers the degree of any logarithmic series by one. We now demonstrate 
via the following string of identities that Q is D-invariant: 

b 

= ELfcl V a -^{ EZ pf\ X )\ 0) Pa- b -M 



b 



E 



a-l 
b 



Lal(^ 0) (^L,K-6-i^) 



(o) J 



E z Q P a a (x). ■ 



We summarize the results of the preceding sections in the following theorem: 
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Theorem 3.6.4 For any graded sequence of formal power series of logarithmic type Pa( x ), the following 
statements are equivalent: 

1. p"(x) is a Roman graded sequence. (Definition 3.2.2) 

2. Pa{x) is the (rfi 1 ) associated graded sequence for some unique delta operator /(D). (Definition 3.3.2j 

3. p"(x) is the (mM 1 ) conjugate graded sequence for some unique delta operator g(D). (Definition 3.5.1J 
4- Pa( x ) is the (nth) fr as j c graded sequence for some unique delta operator /(D). (Definition 3.4. 1) 

5. Pa( x ) is a graded sequence of logarithmic binomial type. (Definition 3.6.2J 

6. Pa( x ) is a graded sequence of formal power series of logarithmic type whose coefficients 

k<n 

satisfy equation (3.4-). 

Furthermore, in this case, p™ (x) is associated with and basic for the same delta operator /(D) (and integer 

n). 

Proof: Condition 1 is equivalent to Condition 6 by Porism 3.2.3, it is equivalent to Condition 3 by 
Proposition 3.5.2, and it is equivalent to Condition 2 by Proposition 3.4.6. Next, Condition 2 is equivalent 
to Condition 4 by Theorem 3.4.2. Finally, Condition 4 is equivalent to Condition 5 by Theorem 3.6.3. 11 



Chapter 4 



Relations Among Roman Graded 
Sequences 

Continuous In this chapter, we focus our attention primarily upon standard operators, and their corre- 
sponding Oth associated and conjugate sequences which we have seen (Propositions 3.3.5C and 3.5.3C) 
are standard sequences. 

Discrete In this chapter, we maintain the same level of generality as before. However, our results are still 
less general than the corresponding results in the continuous case. 



4.1 Transfer Operators 



Let p"(x) be the standard associated graded sequence for the delta operator /(D). In view of Propositions 
3.5.2, and 3.4.6, the standard conjugate graded sequence for the operator /(D) is in general some other 
Roman graded sequence, say g"(x). By Proposition 3.3.1, this graded sequence is in turn the standard 
associated graded sequence for another delta operator .g(D). What is the relationship between /(D) and 
5(D)? And between p" (x) and g(D)7 

We shall prove (Corollary 4.2.9) the remarkable fact that the formal power series /(D) and 5(D) are 
inverses to each other in the sense of LaGrange Inversion [25], and p"(x) is the conjugate graded sequence 
for 5(D). Actually, the results we shall obtain arc more sweeping, and lead to a powerful technique for 
establishing identities among formal power series of logarithmic type. 

We shall occasionally use the boldface p to denote a logarithmic graded sequence p"(x)- Thus, 
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Definition 4.1.1 (Umbral Composition) Let p and q be two graded sequences, and let c a b denote 
the coefficients of the p sequence c a b — [A^ {x)]p™ (x) . We define the umbral composition of the q graded 
sequence with the p graded sequence to be the graded sequence defined by the following convergent summation 

<ta(p)=^2c abPb *(x). 

b 

Proposition 4.1.2 Given two graded sequences p"{x) and q™{x). Their composition p"(q) is a well 
defined graded sequence. 

Proof: We need only consider the continuous case. By the Noetherian condition, their composition is well 
defined. The conditions on its degree and order follow immediately. Continuity follows since the composition 
of two continuous functions is itself continuous. 3 

By [25], D is the two-sided identity for the group of delta operators under (O-)composition, by definition 
the graded sequence of harmonic logarithms is the two-sided identity for the semigroup formed by the 
operation of composition of standard Roman graded sequences. We shall show that this semigroup is actually 
a group, and that the two groups are naturally isomorphic. The isomorphism is given by the function which 
associates each delta operator with its associated graded sequence. The crucial role in obtaining these results 
is played by the notion of a transfer operator which we proceed to define: 

Definition 4.1.3 (Transfer Operator) Let p%(x) be a Roman graded sequence. The transfer operator 
associated with the graded sequence is the continuous linear operator t p : T — > X defined as 

T v K{x)=p a a {x) 

for all a and a. 

Thus, t v <£{x) = g°(p). 

Discrete If p"(x) is the associated graded sequence for the delta operator /(D), we frequently write r/ for 
the transfer operator associated with p"{x). 

Continuous If p"{x) is the n^h associated graded sequence for the delta operator ./(D), we frequently write 
Tf. n for the transfer operator associated with p"(x). We also write r/ for r/ ; o- 

Proposition 4.1.4 The transfer operator is a well defined regular operator. 

Proof: Proposition 4.1.2. 1 
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4.2 Adjoints 



We must now define the adjoint of a regular operator. This definition is equivalent to the usual adjoint 
with respect to the inner product defined in §4.3.1 A. 

Definition 4.2.1 (Adjoint) If 9 is a continuous linear operator onT a for some a ^ (0); then the adjoint 
of 9 is defined to be the linear operator adj(#) on Artinian operators defined by 

([adj(0)/(D)] p(x)) a = (f(V)6p(x)) a (4.1) 

for all formal power series of logarithmic type p(x) e l a , and for all Artinian operators /(D) e A+. 

Discrete In the discrete case, if 9 is a continuous, linear operator onl^ ', then the adjoint of 9 is defined 
to be the linear operator adj(0) : A — > A such that equation (4-1) holds for all formal power series of 
logarithmic type p(x) <G and for all differential operators /(D) e A. 

Proposition 4.2.2 Let 9 be as in Definition 4.2.1. Then adj(#) is a well defined continuous linear 
operator. 

Proof: (Well Defined) Spanning argument (Proposition 2.5.5). 

(Linear) Consider the following string of equalities 

(adj(ai0i + a 2 9 2 )f(B)p(x)) a = (/(D)(a 1 1 + a 2 9 2 )p(x)) a 

= ai (f(D)6 lP (x)) a + a 2 (f (D)9 2 p(x)) a 

= ai (adj(0i)/(D)p(aO> a + a 2 { & &](9 2 )f{T))p{x)) a 

= ((aradj^i) + a 2& d](9 2 ))f{D)p{x)) a . 

(Continuity) Note that the expression (f(D)9p(x)) a = (adj(0)/(D)p(x)) Q , is continuous in /(D) and 
p(x). I 

Let 9 be a continuous linear operator on X or T + which maps each I a to itself. The adjoint of the 
restriction of 9 to I a is denoted by adj(0) a . If 9 is a regular, continuous, linear operator on I + or all of J, 
the adjoint of 9 is the unique operator which coincides with adj(0) a for a ^ (0). 

The adjoint of nonrcgular operators on 1 or 1 + can be similarly defined; however, we omit the discussion, 
since it would be an unnecessary digression. 

Let us consider several important operators, and as an exercise compute their adjoints. 
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Example 4.2.1 Let /(D) be an Artinian operator. Then 



[adj(/(D))] ( 5 (D)) = /(D)g(D). 



In other words, the adjoint of an Artinian operator /(D) is the operator of multiplication by /(D). 

The Roman shift is an operator which is crucial to our work in the remainder of this chapter; moreover, 
it is not a D-invariant operator, so its adjoint is of particular interest. 

Example 4.2.2 We define a regular, linear operator a onl continuous on eachl a by requiring that for 
all a, 



We call the operator a the standard Roman shift. It is not a D-invariant operator. For example, Dcr — ctD = 
1^0. The standard Roman shift has as its adjoint the Pincherle derivative operator adj(cr)(D a ) = aD a_1 
for all integers a, since 



Continuous In the continuous case, a is not continuous over all of l a since the limit of ap{x) as p(x) 
approaches 1/x is not well defined. 

The standard Roman shift is a logarithmic generalization of the operator x of multiplication by x, since 
ax a = x a+1 for a ^ —1. See Theorem 4.3.7 for an amazing property obeyed by the Roman shift. 

The Pincherle derivative has several equivalent definitions. 

Definition 4.2.3 (Pincherle Derivative) Let /(D) be an Artinian operator. We can define its Pincherle 
derivative /'(D) by any of the following equivalent formulations: 

1. The Pincherle derivative is the continuous, linear map 



3. The Pincherle derivative is the adjoint of the standard Roman shift a. In other words, /'(D) = 




(bD^X^x 



)) a = [b-]\(l-S , b )S a+1 , b = (D b a\^x 



)>«■ 



A+ -» A+ 
D a i — ► aD 11 - 1 . 



2. 



/'(D) = /(D)a-<7/(D). 



adj(a)/(D). 



4.2. AD JOINTS 



41 



4- The Pincherle derivative of an Artinian operator is its derivative as an Artinian series. (See the section 
on derivatives in [25], and Theorem 2.4.1 here.) 

Before considering the adjoint of an arbitrary transfer operator. Let us consider the following transfer 
operator which illustrates an interesting quirk of the continuous theory; it clarifies the meaning of the "; n" 
in the definition of the exponentiation of Artinian series, and their composition. 

Example 4.2.3C Let tp n be the transfer operator ttj.„. Thus, i[) n is the continuous linear operator 



Finally, we consider the adjoint of an arbitrary transfer operator. 

Proposition 4.2.4 If t is a transfer operator, then its adjoint adj(r) is an automorphism of A + . 



Proof: We show that adj(r) is an monomorphism and that it acts on delta operators. This implies that 
adj(r) preserves degree, and thus, that adj(r) is an automorphism. 

Let t be associated with the Roman graded sequence p" (x) . 

(Injectivity) Assume that adj(r)/(D) = adj(r)g(D) for some pair of Laurent operators /(D) and (/(D). 



so by the spanning argument (Proposition 2.5.5), we infer /(D) = g(D). 

(Morphism) By Proposition 4.2.2, adj(r) is continuous and linear, so we need only confirm that adj(r) 
preserves multiplication. Let /(D) and g(D) be Artinian operators. Then we have for all a, 




Discrete Moreover, adj(r)^ 



is an automorphism of A. 



Thus, for aRp{x) Gl (1) , 



(f(B)rp(x)) {1) = (g(D)Tp(x)\ 



(i)' 



(adj(r)(/(D) 3 (D))A( 1 )(x)) ( 



(i) 



(f(D)g(D)T\W(x)\. 
(/(D). 9 (D)rf)(x)) (1) 



(i) 




b 




b 



E [I] (ladj(r)/(D)] \£\x)) (1) ([adj(r). 9 (D)] A« b (x)) ( 

b 

([adj(r)/(D)] [adj(r) 3 (D)] X£Hx)) (1) 



(i) 
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Thus, adj(r) [/(D)g(D)] = (adj(r)/(D))(adj(r)fif(D)) by the spanning argument. 



(Degree) We need only consider the continuous 
sequence for the delta operator /(D), then for all a, 



case. 



Suppose p%(x) is the (n* n ) associated graded 



(adj(r)/(D)^Ai 1 )(a ; )) ( 



(i) 



[a] \S ab 



(i) 



and by the spanning argument (Proposition 2.5.5), adj(r)/(D) = D. The conclusion now follows by the 



The most important properties of transfer operators are stated in the following proposition. 

Proposition 4.2.5 1. A transfer operator maps Roman graded sequences to Roman graded sequences. 
2. If t : p"(x) i — ► Qa( x ) * s a continuous linear operator, where the p"(x) and q^{x) are the (rfi 1 and 



3. Continuous Moreover, adj(r)<7(D) a; " = /(D) b;n for all a, and if p"(x) and q"(x) are standard then 
the operator t above is a transfer operator. 

Discrete The operator t above is a transfer operator. 
Proof: We need only consider the continuous case. 

(1) Let t : X"(x) i— > p"(x) be a transfer operator. By Theorem 4.2.4, adj(r) is an isomorphism of A + . 
Let q%(x) be the mth associated graded sequence for the delta operator g(D). Then, 



Expansion Theorem 2.5.4. 



I 




(adj(T)-VD) 6;ro T<£(z)>a = <$(D) 6!m £Cr)>a - l^V-Sab- 

By [25] and Proposition 4.2.4, adj(r)- 1 3 (D) b;m = (ad^r)- 1 g(D)) b '> n for some n. 

Hence, rq^(x) is an associated graded sequence for the delta operator adj(T) _1 <7(D). 



(2) We have the following sequence of equalities: 



(adj(r) 3 (D)rf)(x)) ( 



(i) 



(g(B)rp^(x)) (1) 
(g(D) q W(x)) {1) 

Sa.l 



(f(D)p£Xx)) { 



4.2. AD JOINTS 



43 



Hence, by the spanning argument, adj(r)g(D) = /(D). 

(3) Similar to (2). 

(4) More generally, for an arbitrary Artinian operator ^ b ct,g(D) b ' ,n , we have by Proposition 4.2.4 and 
[25], adi(T)^2 b c b g{T>) b < n = ]T Q,. 9 (D)<>< m for some m. Specialize to the case of £ & c 6 D b = /( _1 )(D) 6; ™ to get 

adj(r) (/(- 1 )( 5 )) 6i " = e 2 ' r * m D 6 ) 

and hence 

((.f^Hg^rKix))* = (e 2 ™ n D b \:(x)) a . 

In other words, T e 27Tlam A"(x) is associated with the delta operator f^ _1 \g). Finally, note that m — when 
all of the sequences in question are standard. 1 

The following results illustrate the applications of the preceding proposition. 

Proposition 4.2.6 Ifp"(x) and q"(x) are the standard associated graded sequences of the delta operators 
g(D) and /(D) respectively, then the composition g(f) (resp. g(f;0)) is the delta operator with standard 
associated graded sequence q"(p)- 

Proof: If t : X"(x) i— ► p"(x) is a transfer operator, then as in the proof of Part (1) of Proposition 
4.2.5, rq"(x) is the standard associated graded sequence for (adj(r)) _1 .g(D). However, rq"(x) = q" (p) by 
Definition 4.1.1. Moreover, Part (2) of Proposition 4.2.5 asserts that adj(r)/(D) = D, so adj(r) _1 D = /(D), 
and, thus, adj(r) _1 g(D) = g(f) (resp. <?(/;0)). The conclusion follows. I 

Corollary 4.2.7 The set of Roman graded sequences is closed under umbral composition. Similarly, 
the set of standard Roman graded sequences is closed under umbral composition. II 

COROLLARY 4.2.8 (Inverses) Letp"(x), andq"(x) be (standard) Roman graded sequences associated with 
the (standard) delta operators /(D), and g(D) respectively. Then the following statements are equivalent: 

1. f(g) = B (resp. f(g;0) = B), 

2. «?(/) = D (resp. f(g;0) = D), 

3. <£(p) = X°(x), and 

I p«(q) = >%(x). I 

Corollary 4.2.9 The standard associated Roman graded sequence for the standard operator /(D) is 
the standard conjugate Roman graded sequence for the delta operator /( -1;0 )(D) , and visa versa. 



Proof: Theorem 3.3.4 and Corollary 4.2.8. 



■ 
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4.3 The Roman Shift 



The following definition generalizes to all Roman graded sequences the notion of the standard Roman 
shift (Example 4.2.2). 



Definition 4.3.1 (Roman Shift) Ifp"(x) is a Roman graded sequence, the roman shift relative to p"(x) 
is the linear operator a p :2^I defined by 



<rM*) = \ n t 1 (4-2) 

{ ifa = -l. 

for all a and a where cr p is continuous over each l a but not continuous simultaneously over all of I. 



As with Transfer operators, if p"(x) is (n-)associated with the delta operator /(D), we also write aj (or 
<7/ ; „) instead of a p . When no graded sequence of delta operator has been specified, we assume a = crr> = o~x, 
that is, that a is the standard Roman shift as previously defined. 



Definition 4.3.2 (Artinian Derivation) Let Q be an operator on the Artinian (resp. Noetherian) algebra 
such that 

Qf(x) a ' n g(x) b ' m - a(Q/(x))(/(x) a - 1 ^. 9 (x)^ m + 6(Q.g(x))/(x) a -" 3 (x) b - 1 ' m . 



Note that an Artinian (resp. Noetherian) derivation is automatically a derivation, and thus is linear. 
Lemma 4.3.3 The derivative is an Artinian (resp. Noetherian) derivation. I 



Proposition 4.3.4 1. A regular, continuous, linear operator 9 defined on I + is a Roman shift if 
and only z/adj(#) is a continuous, everywhere defined (Artinian) derivation of the algebra of Artinian 
operators A + for which adj(0/(D)) = 1 for some delta operator /(D). 

2. A regular, continuous, linear operator 9 defined on the logarithmic algebra I is a Roman shift if and 
only if &d]{9) is a continuous, everywhere defined, (Artinian) derivation of the algebra of differential 
operators A for which adj(#/(D)) = 1 for some delta operator /(D). 



4.3. THE ROMAN SHIFT 



45 



Proof: We need only consider the continuous case. 

(Only if) Let p"(x) be a Roman graded sequence of formal power series of logarithmic type. By 
Proposition 3.5.2, p"(x) is (n— ) associated with a delta operator /(D). Now, we have 

(adj(a / )/(D)^>i 1 )(x)) (i) = (/(D)%p^W) (1) 

= (a + 1) [a] \6 a +i,b 

= H a 1 !(5 a,6-1 

Hence, by the spanning argument, adj(cr/)/(D) &;Tl = &/(D) b_1; ™. Also, adj(<7/)/(D) = 1. By the continuity 
of adj(<7/), the result follows. 

(If) Conversely, suppose adj(#) is a continuous, everywhere defined derivation of A + with adj(0)/(D) = 1. 
Let <t p be the Roman shift associated with p%(x), the ra tri associated graded sequence for /(D). We shall 
show that 8 = a p . 

(f(Dr n e P £\x)) {1) = (adj(0)/(D)^>W(x)) (1) 

= (a/(D) a - 1 ^>«( a; )) (1) 
= alb~}\5a : b-i 

Thus, by the spanning argument, 9 = er p . 1 
Next, we derive the chain rule for Roman shifts. 

Proposition 4.3.5 (Chain Rule) Suppose (jf, and a g are Roman shift operators. Then 

adj^/) = (adj(cr / ) 5 (D))adj(cr 9 ). 

Proof: We need only consider the continuous case. 
For any Artinian operator h(D) = J^b c fc.9(D) b; ™, 

adjfo)fc(D) = ^& C6 .g(D)"- 1 ^adj(a / ) 3 (D) 

b 

= [adj( CT/ ) 5 (D)] [adj(a fl )/i(D)] , 

so 

adj(cr / ) = (adj(cr / )5(D)) adj(<7 g ). ■ 
The following proposition allows us to relate two Roman shift operators. 
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Proposition 4.3.6 Ifo-f, and a g are regular shift operators, then 

o-f = f7 s adj((7 / )fif(D). 

Proof: For any Artinian series h(D) and logarithmic series p(x), we have 

h(D)a f p(x) = adj{a f )h(D)p(x) 

= (adjK)/ i (D))(adj(a / ) 5 (D))p( a; ) 

by the chain rule. This in turn equals h(D)a g (adi(af)g{D))p(x). I 
4.3.1A Orthogonal Functions 

We can now define a symmetric non-degenerate bilinear form associated with any Roman graded sequence 
Pai x )- Consider the principal subsequence p a (x) (Definition 3.1.1). Define an inner product 

(Pa(x)\p b (x)) = S ab [a]l. 

One verifies that relative to this bilinear form, the operator /(D) is adjoint to the Roman shift 07. For 
example, for the harmonic logarithm one has X a (x) = x a for all a, and 

(Dx a \x b ) = (x a \ax b ) 

where a is the standard Roman shift, which restricts to the operator x of multiplication by x. The principal 
sequence p a (x) is the set of eigenfunctions of the operator cr//(D) with eigenvalue n. For example, the x a 
form the set of eigenvalues of the operator xD. 

The bilinear form is not definite, since (l + x~ 2 |l + x~ 2 ) = 0. 

Discrete However, in the discrete case, we can define a Hermitian form 

(p(x)\q(x)) c = - (p(x)\q(x)}, 

so we have 

(wiw B )c = (Wri^xr) 

{— (x n \x n ) for n even, and 
(x n \x n ) for n odd, 

and thus for n negative 

((ix) n \(ix) n ) c = l/(-n- 1)! > 0. 
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Extending by linearity, we obtain a Hermitian inner product which is positive definite, where the 
sequence h n (x) defined as 

f lx n for n > 0, and 



h n (x) 



(ix) n for n < 



is a complete orthogonal sequence in the Hilbcrt space obtained by completion relative to this Hermitian 
inner product. We can obtain this sort of sequence using the Knuth coefficients [19] with e = — i in 
place of the Roman coefficients. In this manner, 



n\ if n > 0, and 

(-l) n i/(-n- 1)! ifn<0. 



One can therefore develop a spectral theory of the operator xD in this Hilbert space (rather than with 
the smaller one including only positive powers of x, as is done classically). 



4.3.2A Rodrigues' Formula 



Note the following amazing fact. 



Theorem 4.3.7 Let p(x) be a finite linear combination of harmonic logarithms X^ =1 CjA"' ' (x). Then 
for any complex number z, we have the following formal identity 

(e- za T>e za )p{x) = (D - zl)p(x). 



By e acr De ao p{x) we merely mean to indicate the expression 

x x ^ (-l)V+ fe . ^ k 
j>0fe>0 J ' 

which we assert does in fact converge to the indicated value. The "operator" e za can not be extended to all 
of X, since the series X)fe>o zk(jk A- 

The corresponding classical identity 

e" x De x = D + I 

is associated with the classical identity 

Dx - xD = I 

which corresponds to the logarithmic identity 

Dcr - crD = D' = I. 
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Proof: By linearity, it suffices to consider the case p(x) = \"{x). There are three cases to consider. 
(When a is not a nonpositive integer) Classical proof can be applied mutatis mutandis. 
(When a = 0) We calculate: 

e~ za De zcr Xq(x) = ^Et^W 
j>o fc>o ■'■ " 

j>0fe>0 J ' 



k>0 J ' \fc=0 u/ , 



= A^-zA^x) 
= (D-zI)A«(x) 



(When a is a negative integer) We similarly compute 

l) j a 



e— De-A^(x) = EEt^^'^W 



j>0 fc>0 
fe>0j=0"'' v 

V^ 1 V- (-f) J z fc (fl + fc-j) 

fe=o i=o ■'^ 
= aX^x) + zX^x) 
= (D-zI)X:(x). ■ 



4.4 Explicit Formulas for Roman Graded Sequences 

We are now ready to derive the following recurrence formula for Roman graded sequences: 

Theorem 4.4.1 (Recurrence Formula) Ifp"(x) is an associated graded sequence of delta operator /(D), 
then for all a ^ —1 and for all a 

p a a+1 (x)=a(f'(D))-' P :(x) 

where a is the standard Roman shift. 
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Proof: Let a^-1, then by Proposition 4.3.6. 

Pa+l( X ) = °fP2( X ) 

= ( T(adj( ( T).f(D))- 1 p^(^) 
= af'(D)~ 1 p"(x). I 



Next, we give an explicit formula for the associated graded sequence of a delta operator in terms of the 
residual series of the graded sequence of harmonic logarithms. 

Proposition 4.4.2 Ifp^(x) is the (vM 1 ) associated graded sequence of formal power series of logarithmic 
type for the delta operator /(D), then for all a and all a ^ (0): 

f Discrete p°{x) = [a] !/'(D)/(D)~ 1 ~°A" 1 (a;) 1 
1 Continuous p a a {x) = [a] If (D)/(D)- 1 -°' n A° 1 (a;) J 



Proof: We need consider only the continuous case. Let q"{x) be the graded sequence defined by: 

qZ(x)=la]\f'(D)f(V)- 1 - a \1 1 (x). 

It suffices to verify that q% (x) is the n^h basic graded sequence for /(D). q%(x) is indeed is a graded sequence 
since deg (/'(D)/(D) _1 — ; ") = — 1 — a and all the operators in question are continuous. Then, if a ^ 0, we 
have 

<£(*)><, = (H!/'(D)/(D)- 1 -^A« 1 (x)) Q 

= ^((/(Dr^'A-^ 

Now, ((/(D)-)' A^(x)) Q = [D- 1 ] (/(D)-)'. 

The crucial observation in any theorem related to Lagrange inversion is that for all Artinian operators 
(/(D), the coefficient [D _1 ]<7'(D) is zero. Hence, {qa(x)) a = 0. Thus, Property 2 of Definition 3.4.1 is 
satisfied. 

Now, consider the case, a = 0, and (0). Then we have: 

(9oW) Q = (/'(D)/(D)- 1 A« 1 (x)) Q . 

Say that c b = [D b /(D)]. The coefficient [D" 1 ]/^" 1 = c^ 1 , and the coefficient [D°]/'(D) = Cl . Neither 
operator has any terms of lower degree. Hence, the coefficient [D -1 ] (/'(D)/(D) _1 ) = 1, and there are no 
terms of lower degree. Hence, (q§ (x)) a = 1. Thus, Property f of Definition 3.4.1 is satisfied. 
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Finally, 

/(D)*i"(£(aO = [a]\f(D)f(D)- a X1 1 (x) 

so Property 3 is also satisfied. Therefore, p%(x) — q%(x). B 

Proposition 4.4.2 is unusual in that it does not readily generalize to 1^ because for a = (0), A" 1 (x) 
equals 0. Moreover, for a > 0, /(D)~ 1_a (resp. /(D) _1_a ; n) has negative degree, and thus is not a member 
of A. 

Nonetheless, its consequence — the transfer formula — still holds for T^: 



Theorem 4.4.3 (Transfer Formula) Let /(D) be a differential operator of degree one, and let p"(x) be 
its (vft 1 ) associated graded sequence. Then for all a and a, 

Discrete p°(x) = /'(D) ( + X"(x) 



Continuous p a a {x) = f'(p)[JL\ ' X"(x) 



/(D) J 

— Y 

In particular, its residual series is given by 

pV(x) = f'(D)x- 1 . 

Proof: (When a ^ (0)) The conclusion is immediate from the preceding proposition. 

(When a = (0)) We need only consider the continuous case. By regularity, we have the following string 
of equalities for all a: 

= ^(o),(i)/'(D) [j^ j A«(x) 
= ^( D )(/(5)J E {0U1) X^(x) 

a+l;n 

I 



f (D) (^) »?»«. 



The following variants of the transfer formula arc often useful: 
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Corollary 4.4.4 Letp"(x) and q"(x) be the (rfi 1 and mM 1 ) associated sequences for the delta operators 
/(D) and g(D) respectively, then for all a and a, 

| Discre «e rfW . 1 

PROPOSITION 4.4.5 In the notation of Proposition 4.4.2, 

Discrete p<£(x) = fl (D)-A«(a;) - (a(D)- a )' X^x) 1 
Continuous p°(x) - g(D)~ a; "A"(ir) - (,g(D) _a;I1 )' A"_ x (x) J 

/or a ^ 0, w/iere 5 (D) = /(D)/D. 

Proof: We need only consider the continuous case. By Proposition 4.4.2, 

P :(x) = f'(B)g(B)- 1 - a -' n X:(x). 

However, 

/'(D)g(D)- 1 - a = (D.g(D))' 5 (D)- a - 1 ^ 

= (D'. 9 (D)+D 3 '(D))g(D)- a - 1 ^ 
= .g(D)- a + D. 9 '(D).o(D)- a - 1 ^ 
= . (D)- a + (.g(D)-^)'D/a 

so that 

f'(D)g(D)- 1 - a -' n X a a (x) = g(D)-^X a a (x) + (ff(D)"^)' A^(x). ■ 
Note that Proposition 4.4.5 does not in general hold for a = 0. 

Corollary 4.4.6 In the notation of Proposition 4.4.5, 

Discrete p a a {x) = ag(D)~ a X"_ 1 (x) j 
Continuous p°{x) = ag(D)^ a:n X"_ 1 (x) J 

for a^0,l, where a is the standard Roman shift. 
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Proof: We need only consider the continuous case. By Proposition 4.4.5, 

= g(D)- a -' n K(x) g^y^aX^x) + <jg(T))~ a ' n X™_ 1 {x). I 

We conclude with an important remark about graded sequences of formal power series of logarithmic 
type. Let be a Roman graded sequence with coefficients ci"f = [A^ (x)]p% (x) . By regularity, 

1. For a, [3 + (0) and for all a and b, d%£ = dff , 

2. For a and b not a negative integer, df,j^ = d^°\ and 

3. For a ^ [3, and for all a and b, = 0. 

In view of this, we see that in computations with Roman graded sequences it suffices for most purposes to 
compute in the subspace In other words, even in computations with polynomials it is preferable to deal 
with logarithms first! A quick survey of the examples chapter 5 demonstrates the utility of as compared 
the algebra of polynomials. 



4.5 Composition of Formal Series 



Let /(D) be a delta operator. All of the formula for composition and inversions of series (in particular, 
the various versions of the LaGrange inversion formula) are consequences of the following theorem. See also 
[25] and [22] for other results concerning the composition of series. 



Proposition 4.5.1 1. // /(D) is the delta operator with (n™) associated graded sequence p^(x), then 
for every Artinian operator g(D) we have the following convergent series: 

' Discrete - VJ«W, D » 1 

Continuous «(/<-'«) = ^MEMilko' ' 

b 

where a / (0). 

2. equation (4-3) holds for a = (0) as well whenever g(D) is a differential operator. 
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Proof: We need only consider the continuous case. By the expansion theorem (Theorem 3.3.3), 

g (D)=E to( y >a /(D)^. 

Hence, substituting /( -1; ™)(D) for D (by the characterization of composition in [25]), 

= ^ (.9(DK(x)) a /(/( _ 1;Tt)))fc; „ 



Theorem 4.5.1 has many versions and many corollaries. We first deduce from the Expansion Theorem 
the following convergent expansion 



N! 

Continuous = g {f ^) = ^M!3^D t 



N! 



L6]! 



L6]! 



Hence, 

Discrete (ff(D)p°(x)) Q - (g(f^)X^(x) ] 
Continuous {g(T>)p%(x)) a - ($(/ ( ~ lin) )Ai (z); 
. An application of the Transfer Formula (Theorem 4.4.3) gives (in, for example, the continuous case) 

Lal!( 3 (D)/'(D)/(D)- 1 -^Ai 1 )(x)) (1) = (g(D)pW (xj) (1) 
By the spanning argument (Proposition 2.5.5), we have the following corollary. 



Corollary 4.5.2 Let /(D) be a delta operator, and let g(D) be any Artinian operator. Then we have 
the following convergent sum: 

Discrete g(f^) = ]T (. 9 (D)/'(D)/(D)- 1 - fe A«( a; )) (1) D fc 

k 

Continuous 5 (/ ( " 1;n) ) = E (s(D)/'(D)/(D)- 1 -°" l A^(a;)) (1) D <l . I 



By taking g(D) = D a , we obtain powers of /("^(D) (resp. /("^(D)). 
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COROLLARY 4.5.3 ///(D) is a delta operator with (rfi 1 ) associated graded sequence p"(x), then we have 
the following convergent expansion 



Discrete / (_1) (D)° 
Continuous /( _1 ; n )(D)° = ^ 



Chapter 5 

Examples 



5.1 Roman Graded Sequences 



We prefix some general considerations about the computation of the Roman graded sequence p% (x) (n- 
)associated with a delta operator /(D). The crucial step is the computation of the residual series p^\(x). 
This is given by the simple formula 

p W 1 (x)=f'(D)(l/x). 

Once the residual series is known, any of the scries Pa (x) can be obtained from the residual scries by 
applying a suitable power of /(D); that is, by the Transfer Formula (Theorem 4.4.3), 

/ fmn - fe " 1 

Discrete p%(x) = /'(D) 



Continuous p%(x) 
when a is a negative integer, and 

Discrete p%(x) = /'(D) 

Continuous p°(x) = /'(D) 
otherwise. 



/'(D) 



X 

-a—l;n 



-k-l 



x k |^log x-l 

-a— l;n / 

x a I log X — 



1 

~ 2 ~ " 
L-a]! 
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If we denote the coefficients of p a L> {x) by c a b = [A^ {x)]p a L> (x), then it follows from regularity that 

b 

for all a. 

Note that for a = (0) we obtain 



Discrete j4 (x) = ^c^-x- 7 

Continuous p a ° a \x) = ^ c a 

&eR-N- 



In the discrete case, this is the sequence of polynomials of binomial type associated with the delta operator 
/(D); thus we see that even in the case of polynomials, it may be speedier to compute via the logarithmic 
graded sequence. 

Table 5.1: Examples of Roman Graded Sequences 



Delta 
Operator 



Associated 
Graded Sequence 
TO 

(*)n 

c;w 



Inverse 
Operator 



Conjugate 
Graded Sequence 



Qn( x 

Ki x ) 



D 

A = E - I 
V = I-fi- 1 
A a = DP 
E Z (E - 1) 

K = D(D-I) 



D 

log(I + D) 
-log(I-D) 



K 



D/(D-I) 



5.1.1 Lower Factorial 



Other than the Harmonic logarithms (which we have already discussed at great length), our first example 
of a Roman graded sequence is the logarithmic lower factorial graded sequence. 



Definition 5.1.1 (Forward Difference Operator) Define the forward difference operator A = E — I = 
e D — I. Let (x)" denote its standard associated graded sequence; it is called the logarithmic lower factorial 
graded sequence. Let </>"(x) denote its standard conjugate graded sequence; it is called the logarithmic 
exponential graded sequence. 
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Table 5.2: Lower Factorials (x)" 



(4 o) = 


= x(x — 1) 


(X)' 1 ' = 


= \i 1} (x) - Aj 1 ' (x) + B 3 , 3 /3x - B 4 , 3 /12x 2 + ■ • ■ 


(x)< > = 


- X 


(*) 1} = 
{*h - 

(*)<_!> = 


xlog(x) - a: + B 2 , 2 /2x - B 3 , 2 /6x 2 H 


(4 o) = 


1 


= log(x + 1) + Bi/(1 + x) - B 2 /2(l + x) 2 + B 3 /3(l + x) 3 

= l/(x + l) 

= l/(x + l)(x + 2) 



Now, A' = £7, so by Theorem 4.4.3, we calculate the residual series 

(D - BV"1 - _JL_ 



{ x y_{ = E X -' = 



x + 1 



In general, 

PROPOSITION 5.1.2 For n a negative integer, 



1 



(x + 1) • • • (x — n) 



{x) n - 



(5.1) 



Proof: We proceed by induction. The case n = 1 amounts to equation (5.1) which we verified above. 
Now, suppose proposition holds for n then 



( x )-n-l 



[-ny'Aixfl 
If 1 



1 



n \{x + !)■■■ {x + n) (x + 2) • • • (x + n + 1) 
1 



(ar + l)---(x + n+l)' 
Similarly, we have the continuous analog of [39, p. 133]: 

PROPOSITION 5.1.3 For a not a negative integer, (x)^ = (x) a . 

Proof: Proposition 3.3.1 1 
Thus, 



Corollary 5.1.4 For all a, (x) a = (x) a . 
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Now, we determine the series (x)^ . By Theorem 4.4.3, we have 

{x)P = E^logx, 

or equivalcntly 



I 



(t)^dt = log(x + l). 



Moreover, since by the Euler-MacLaurin formula 



^=E B * (5-2) 



fe>0 

we have 

(-)o 1} = Ef^k^ + i) 

k>o K - 

where the are the Bernoulli numbers. Hence, we discover that (x)^ = ip(x + 1) coincides with the 
classical ip- function (the logarithmic derivative of the gamma function) introduced by Gauss. Similarly, one 
finds that (x)^ and (x)^ coincide with the digamma and trigamma functions. 



The classical expansion 

Discrete (D/A)™ = ^B kn T) k /k\ 

fe>0 

Continuous (D/A) Q;0 = ^B ka V k /k\ 



fe>0 



defines the graded sequence B a b of Bernoulli numbers of order b and degree a. In terms of these higher 
order Bernoulli numbers, we obtain for all a: 

(x): = (D/A) 1+a;0 EX"(x) 

= Y,B k , a+1 n k /k\\%(x + l) 



fe>0 



5> M+1 a k ]x:_ k (x+i) 



k>0 



(2) 

Discrete We now give an explicit calculation of for n a negative integer. To begin with, by Theorem 

4.4.3, we can calculate the residual series of order (2). 

{x) { *\ = E(2x~ 1 logx) 
21og(a; + 1) 
x + 1 
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Continuing in this way, 



= 2 



(x)™-(x + l)™ 

' logjx + 1) logjx + 2) 
x + 1 x + 2 

log(x + l)- (x + l)log(f±f 



(x + l)(x + 2) 



and 



log(x + 1) - (x + 1) log (§±f ) log(x + 2) - (x + 2) log 



(x + l)(x + 2) 
21og(x + l) lo s(f$i) 



(x + 2)(x + 3) 



(x + l)(x + 2)(x + 3) (x + 3) 
In general for n a positive integer, 



21og(jg 
(x + 1) • • • (x + n) ' (n-l)!(x + n)' 



Similarly, we calculate (x)! ' 1 ^. The residual series of order (0,1) is: 

= S(l/xlogx) 

= l/(x + l)log(x + l). 



Next, 



(x + 1) log(f±f + log(x + 2) 



(x + l)(x + 2) log(x + 1) log(x + 2) ' 



and so on. 



The identity 



(x + ay 



U)_ E 

fe>0 
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gives a classical identity satisfied by the -0-function, that is: 

(-l) fc+1 a(a-l)---(a-fc + l) 
, >0 k(x + l)(x + 2)---(x + k) 



i>{x + a+l) = tjj(x + 1) + 



Similar identities can be obtained for the digamma and trigamma functions. 



The logarithmic Taylor's theorem (Theorem 3.3.4) gives the following generalization of Newton's expan- 
sion: 



Proposition 5.1.5 Every logarithmic series p(x) can be uniquely expanded as a convergent series 

p(x)=J2d:(x):/la}\ 

a, a 

where the coefficients are given by 

d: = (A a -<° P (x)) a - ■ 



For example, 



- = J2l/k\(x + 1) • • • {x + k + 1). 



(5.5) 



fe>0 



We digress to indicate the meaning of such equalities. Formally, we merely mean that when 
both sides of the equality are expanded in terms of harmonic logarithms A"(x) the resulting 
coefficients are identical. However, because of such results as Theorem 3.1.3, we are allowed to 
make computations in the real or complex numbers, and thus obtain asymptotic expansions. In 
equation (5.5), the right side and left side are both approximately 0.01 for x — 100; in fact, the 
error is about 0.0000015 when you compute 20 or more terms of the summation. Similarly for 
x = 69, when one computes the first 14 terms of the summation, one finds that the left side is 
0.14488 and the right side is 0.14493. End of Digression. 



Open Problem 5.1.6 How are the Stirling numbers s(a, b) which are the coefficients of the formal power 
series (y) a related to the coefficients of the logarithmic series (x)% ? 
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Table 5.3: Upper Factorials (x}° 



(x) { 2 r) = A< X) (x) + A< X) (x) - B 3 , 3 /3x - B 4 , 3 /12x 2 

(a;)* 1 ' = -xlog(x) + x - Bi^l2x - B 3]2 /6x 2 - ■ ■ ■ 

(xy 1} = log(x-l) - Bi/(x- 1) + B 2 /2(x - l) 2 - B 3 /3(x - l) 3 + ■• ■ 

(x>^ = V(x-l) 

(x> ( _^ = l/(x-l)(x-2) 



5.1.2 Upper Factorial 



Definition 5.1.7 (Backward Difference Operator) Define the backward difference operator V = I — 
E^ 1 =1 — e~ D . Let (x)" denote its standard associated graded sequence; it is called the logarithmic upper 
factorial graded sequence. 

As before, the residual series is given by 
and for n a negative integer 

(x)W=l/(x-l)...(x + n). 
Similarly, for a not a negative integer, we have: 

(x)^=T(x + a)/T(x). 

Thus, for all a, 

(£) a = T{x + a)/T{x) 

For a = 0we have: 

Proposition 5.1.8 (x)^ = log{x - 1) - B x jx - 1 + B 2 /2(x - l) 2 ) - B 3 /3(x - l) 3 + • • • . 
Proof: By Theorem 4.4.3, we have 

(x)^ =E-^\ogx. 
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From the Euler-MacLaurin formula, namely from 

D -D 



V ~ e -D _ i 

= 

= E§(-d)* 

fe>0 

= Em)^- 



k>0 



we infer 



tf> = ^(-l) fe |fD fc log(x + l) 



fc>0 



k) 



\og(x - 1) + 



B, 



x-l 2{x-l) 2 3(a;-l) : 
We have, in terms of the Bernoulli numbers of higher order: 



+ 



Hence, for all a, 



a A) 



D 



k>0 



l+a;0 



= £(-l) fc B fc , 0+1 D fc A2(s-l)/fc! 

fc>0 

= ^(-l) fe S fe , a+1 ^jA^ fe (x-l). 



k>0 

(2) 

Discrete As with the upper factorial graded sequence, we compute for n a negative integer as follows 

starting with the residual series of order (2). 

(x)^ = E- 1 fa-Hog x) 
21og(x- 1) 



x- 1 

v<*><3 

log(x — 2) log(a; — 1) 



x-2 x-l 
log(x-2)-(x-2)k>g(f5±) 
(x-l)(x-2) 
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and in general 

(2) 21og(x-n) 21 °g(j^Q 

[X) - n (x - 1) • • • (x - n) (n - l)\(x - n) ' 



5.1.3 Abel 



The logarithmic extension of the Abel polynomials turns out to be surprisingly pleasing. 

Definition 5.1.9 (Abel Operator) Define the Abel operator A z = DE Z . Its standard associated graded 
sequence A^(x) is called the logarithmic Abel graded sequence. Its standard conjugate graded sequence 
is called the logarithmic inverse- Abel graded sequence. 



Table 5.4: Abel Graded Sequence A%(x) 



4°>(x) 


— x(x — 2z) 


K(x) 
AY\x) 


= crX^Ux - 2z) 


KHz) 


— X 


— (7 log(x — z) 


= 1 


K\x) 
A h \{x) 
A (1 \(x) 


= logO) - zjx 
— x(x + z)~ 2 
= x(x + 2z)~ 3 



By Corollary 4.4.6, for a ^ 0, 1, we obtain 



a-1 



fc>0 

E 

k>0 



k 

a- 1 
k 



{-azfK^ix) 
{-az) k \ a a _ k {x). 



In particular, when a is not a negative integer, we obtain a simple generalization of the classical Abel 
polynomials. 

f Discrete A^\x) = z(i-nz)"" 1 ^ 



1 Continuous A^\x) = i(i-az) a " 1;0 , 
and for a a negative integer, we have 

A { a\x) =x{x-az) a -\ 



(5.6) 



(5.7) 
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Thus, its residual scries is 



(.x + z) 2 ' 
and its principal sequence is 

Discrete A(x) = x/(x - az) a '° 
Continuous A{x) = x/(x — az) a . 

A*q\x) can be calculated via Theorem 4.4.3. It is expressed very simply. 

A { o\x) = A' Z E- Z log x 



= (I + zD)logx 

= logx + z/x, (5. 



so 

A a =\ a {x)+z\ a _ 1 {x). 



From the logarithmic binomial identity, we have the sum 

A i i\x + b) = Y J 



Af[b)A^x) 



k>0 

over all nonncgative integers k. Thus, we infer the remarkable identity 



r+log (x + b = -+log.x + — '— — v 5.9 

x + b by ' x ^ k(x + ak) k+1 



k>l 

For example, we can substitute here the values a = 1, b = 2, and x = 5. If we compute the first 12 terms of 
the series, the left hand and right hand sides of equation (5.9) are both approximately 2.0887673. 

In general, by Theorem 4.4.3 

A«{x) = E- az (l + zD)\ a a (x) 

= \*(x - az) + z [a] \a-i(x- az). 

For example, 

A^ (x) = x log(x — z) + z — x. 

Again, by Theorem 3.3.4, every formal power series of logarithmic type can be expanded in terms of Abel 
series 



n. rv <- ' 
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where 



For example, (See equation (5.13)) 

log a; = ^2(ka) 

k<0 

- aWi 



-k 



= A^'ix) + zA y '[{x) - 2z 2 A [ ^(x) + 9z 3 A W 3 (x) 

fe+i 

= logx V 



fc>o 



x + kz 



z k x. 



That is, 



2-^ \x + ka) 

k>o v 7 



(5.10) 



5.1.4 Gould 



Definition 5.1.10 (Logarithmic Gould Graded Sequence) We define the logarithmic Gould graded 
sequence G%(x) to be the standard graded sequence associated with the delta operator E Z A = E z+1 — E z . 



Table 5.5: Gould Graded Sequence G"(x) 



g' 0) (x) 
G 0) (x) 
G< 0) (x) 


— x(x — 2a — 1) 

— X 

= 1 




G^x) 


- (x> 2 a(x)j +2^ fc > 2 L „_l-|!( x _l)... (x _„ + fc) 




= xiog(x)-x + ^; j , S2 


L<2) i 2 -| i_* . ((!)+-(*- i )(2=i)) ( - 1 >* 


G™{x) 


= ^x-t+Hf-Vt 1 * #(* ^ 2 ^/b!fx 2 ++l)-(x + 2)(x + 3) + • • ■ 


G ( _!>(x) 


— cc/(a; + a) (a: + a + 1) 




G ( _!>(x) 


= + 2a)(x + 2a + l)(x + 2a + 2) 
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The Pinchcrle derivative of E Z A is (z + 1)E Z+1 — aE z , so the residual series is given by 

G ( ±l(x) = ((z + l)E z+1 - aE z ) (- 



z + 1 



X + Z + 1 X + z 
X 

(x + z)(x + z + 1)' 
Since Roman graded sequences are basic, we have 

= E 



x + l 



(x + z)(x + z + l) (x + z + l)(x + z + 2) 

x 



(x + 2z){x + 2z+ l)(x + 2z + 2) 

= xT2-z {x + 2z) -*- 
Similarly, by induction we have for n positive 

G^l(x) = x(x + nz- l)_„_i, 

and, by induction for n nonnegative 

G<£\x)=x(x-nz-l) n - 1 . 

In general for all a, 

Ga(x) = x(x - az - l) -i. 
See §5.2.6 for an explicit computation of G"(x). 

5.1.5 Laguerre 

Our final example of a Roman graded sequence is the logarithmic Laguerre graded sequence. 

Definition 5.1.11 (Laguerre Operator) Define the Laguerre operator to be 

K = D/(D-I). 

Define the logarithmic Laguerre graded sequence to be its (0^) associated graded sequence denoted L"(x). 
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Table 5.6: Laguerre Graded Sequence £"(x) 



L<°>(*) = 
Lf\x) = 

4 0) (-) = 



x* - 2x 
— x 
1 



Lf\x) 

L^l(x) 
r(i) 



= \ ( ^(x)-2\^(x) 

= log(x) - ^ feS0 (-l) fe fe!x- 



L { l\{x) = -J2 k>2 (-l) k (k-Wx- k /2 



Continuous Note that the Laguerre operator K and the Laguerre graded sequence L" (x) are not standard, 
since the leading of X is -1. However, —K and K(—D; 0) are both standard operators we discuss later 
(§5.2.1 and §5.2.2), and if r"(x) and p„ (x) are their standard associated sequences, then 



L«(x)=V-i/ 2 C(x) 
where ij) a is as defined in Example 4.2.3C, and 

l«(x) = (-irx(x). 



(5.11) 



(5.12) 



Now, the Pincherle derivative of the Laguerre operator is given by K' = — (D — I) 2 , so by Theorem 4.4.3, 

L a a (x) = -(D-ir 2 (D-ir +1 ;°A«(x) 

= -(D-ir^x). 



Hence, for all a, 



Thus, 



k>0 ^ ' 

fe>0 



= e 



= e 



L?(x) = -A?(x). 
Similarly, we have derived the following expansion 



hf{x) = log(x) + -- ^ + ^- ^ + . 

X X Z X 4 X 4 



112 6 
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Discrete Thus, L"(x) does not contain any terms of negative degree when n is a positive integer. This 
is true only for Roman sequences associated with D/(aD + 61) where a and b are a nonzero complex 
numbers. Note that the series of degree contains no negative terms only if the operator in question 
is aD where a is a nonzero complex number. 

Continuous In the continuous case, L" contains terms of negative degree for all a. The only Roman se- 
quences which do not always contain terms of negative degree are those associated with delta operators 
of the form zTJ where z is a nonzero complex number. 

From Theorem 4.3.7, we derive a logarithmic extension of the classical Rodrigues' formula for Laguerre 
polynomials. 

LZ(x) = -e°D a - 1 e-°XZ{x). 



5.2 Connection Constants 

Given two graded sequences p%(x) and q% (x); we would like to express one in terms of the other. 

b 

The coefficients d a b are called the connection constants from the graded sequence (x) to the graded sequence 
Pa(x), and are denoted {q£ (x)]p% (x) . 

If Pa{x) and q"(x) are the standard Roman graded sequences associated with the delta operators /(D) 
and g(D) respectively, then by Proposition 4.2.6 

b 

is the standard Roman graded sequence associated with <7(/( _1 );0). Thus, to determine the connection 
constants it suffices merely to calculate r° (x) . This easy device for the computation of connection constants 
is the most effective application of the present theory. 

5.2.1 Upper Factorial to Lower Factorial 

To express (x)" in terms of (x)™ we first calculate 



e -log(I-D) _j 
I 



I-D 

D 
I-D 
= -KID) 



-I 
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where K(D) is the Laguerre operator (Definition 5.1.11). Thus, r"(x) is the standard graded sequence 
related to the logarithmic Laguerre graded sequence defined by equation (5.11), 

where ^1/2 is as defined in Example 4.2.3C 

w\iiK\ x ) = e X a( x )- 
In a sense, V1/2 1S the logarithmic analog of substitution of —x for x. 
We now apply the results from §5.1.5. Thus, for all a and a 



Table 5.7: Lower Factorial in Terms of Upper Factorial 




fe>0 



k J[a-k~\\ 



From Theorem 4.3.7, we have the identity 

f Discrete (x)° = -e <TV V° _1 e- <rv <a;)°. 
\ Continuous (x)° = -e CTv V a - 1;0 e" CTV (x)°. 



5.2.2 Lower Factorial to Upper Factorial 



Conversely, to compute (x)" in terms of (x)", we need the delta operator 

V(A(- 1 '°)) = 5 5_=fT(-D). 

Hence, the connection constants from (x)" to (x)" are given by the coefficients of e~ ma L"(x) (equation 
(5.12)). Hence, for all a and a, 

ws-B-tff 1 ; 1 )^*-.!*). 



fe>0 
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Table 5.8: Upper Factorial in Terms of Lower Factorial 



(x)^ = (x)^ (x) - 2(x)< 1 » 

id 1 ' = (4 1> -E,> (- 1 ) fcfe! ^^- 1 

< X >-2 = E fc > 2 ^ 1 ) fc ( fc - 1 ) fc! ^)^/2 



5.2.3 Laguerre to Harmonic 

If p"{x) is the standard Roman graded sequence associated with the delta operator /(D), then finding the 
connection constants from p% (x) to A" (a;) is tantamount to finding the standard graded sequence associated 
with /( _1; °)(D); that is, the standard conjugate graded sequence for /(D). 

Note that since x/(x—l) is the 0-compositional inverse of itself, the logarithmic Laguerre graded sequence 
is self-conjugate, so we have for all a and a 

fe>0 v 7 L 

and from Theorem 4.3.7, 

K(x) = L°(L) = -e^K-We-'LLZix). 

5.2.4 Lower Factorial to Harmonic 

Similarly, the connection constants from (a;)" to A"(x) are given by <j>"(x) — the logarithmic exponential 
graded sequence — which we now compute. 

The relevant delta operator is 

log(I + D) =^(-l) fe+1 D fe /fc. 

fc>0 

Thus, by Theorem 4.4.3, we can calculate the residual series. 

Ci(z) = (iTd) *-!(*) 

= ^(-l^A^) 

fe>0 

fe>0 L 
fc>0 



<x>°> = (x)°> 
<x>°> = (x)°> 

<4 o) = (x) o) 



2(x)i' 



(0) 
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Finally, by Theorems 4.4.1 and 4.3.7, for a ^ —1, we obtain the recursion formula 

C+i(z) = *(I + D)CW - <7e-De CT CW- 



5.2.5 Abel to Harmonic 



The compositional inverse of the Abel operator is not easily calculated. Nevertheless, we may calculate 
the logarithmic inverse-Abel graded sequence, using the theory of conjugate graded sequences. 



Thus, by Definition 3.5.1 



^ (V b E b *\<* a {x)) a 

fi( X )=2^ p \ 



(x). 



Now, 



Thus, 



B b E bz X2(x) = D b J2 



j>0 



(bzy\«_Ax) 



j>a ■' L J 1 



fc>0 

= £((a-fc)z) fc r]A^ fc (x). 

fe>0 



Hence, the remarkable identity 



fc>0 



(5.13) 
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5.2.6 Upper Factorial to Gould 

Assume that Gould parameter z is a real number which we denote by t. The relevant operator is /(D) = 
— D(I — D)~ t;0 whose associated graded sequence is: 



-a-l;0 

KM 



r») - /'<D)(ffi>y 

= e 2ma {{t- l)D + I)(I-D) at - 1;0 A^(x). (5.14) 



As expected, for t = 1 we rederive the Laguerre graded sequence L"(x), and for t = we get a variant on 
the harmonic graded sequence ^1/2 A" (x). 

For a 7^ 0, 1, instead of equation (5.14), we may use Corollary 4.4.6. 



fe>0 

Thus, for a 7^ 0, 1, we obtain the following remarkable identity relating the Gould graded sequence to the 
lower factorial graded sequence. 

fe>0 

For a = 0, equation (5.14) reduces to 

r?(z) = ((i-lp + IXI-D)" 1 ^*) 

k>0 

Thus, the Gould series of degree zero are given by the elegant expression 

G«(x) = (x)«+tJ2(-l) k+1 (k-iy.(x) a _ k . 

fc>0 

Thus, 



t t t 

+ x + l ~ (x + l)(x + 2) + 2!(x + l)(x + 2)(x + 3) + ' 
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Finally, from Proposition 4.4.5, wc obtain 

r?(aO = -(I-D^Af^-ta-D)*- 1 ^^) 



so 



a ?w + e (Q + *(i _ i)) l 1 - fc i ! 

Af(*)+E((3+*Ql 1 i)) (-!)*(* - 2 ) !A ?-*( a )> 



+ e (Q) + *Q _ 1)) (-!) fc ( fc - 2 ) ! <*>f-* 

(-i) fc ((D+^a- l 1 ))(^-2)! 



Gf(x) 

G^x) = x\og(x) -x + J2B ( ^ I x'- k + J2 



k>2 



k>2 



{x-l)---(x-k + l) 



5.2.7 Lower Factorial to Gould 



Similarly, here we are interested in the operator /(D) = D(I + D)* where t is real. For t = 0, we rederive the 
harmonic logarithms X"(x), and for t = — 1, we get a variant of the Laguerre graded sequence (— l) ma L"(x). 
(See equation (5.12) and §5.2.2.) 



The standard associated graded sequence for /(D) is 



For a ^ 0, 1, 



Thus, for a ^ 0,1, 



-a-l;0 



r a a {x) = /'(D) (^j XZ{x) 

= ((t + l)D + I)(I + D)- at - 1;0 A«(x). 



r»{x) = ail + Dr^X^ix) 



fc>0 
k^a 



k>0 
k^n 
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Now, for a = 0, 

rg(x) = ((t + lp + IXI + D)-^) 

= (1 + tD - tD 2 + tD 3 )A£(x) 

- A«(x)+t]T(fc-l)!A« fc (x), 

fc>0 

so that again 

G%(x) = (x)5(x)+t'E(k-l)\(x)1 k . 

k>0 

For n = 1 , 

rf(x) = (I + D)- t A?(a;)-t(I + D)-*- 1 A^(x) 

= w + E (( fc *) + (- 1 ) fcfc!A -i-*( a; )' 

hence another remarkable expansion for a residual series: 

c°-M=i:(Q +'(1:1)) 

As noted in [20], nearly any sequence may be used as a "factorial" on which to base on umbral calculus. 
In particular, we could have chosen the Gaussian coefficients as our factorials. If so, then would have derived 
the g-analog of our theory. 

Open Problem 5.2.1 What are the q-analogs of these examples? 

We hope the preceding examples display the utility of the theory of formal power series of logarithmic 
type. 
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